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PHYSICAL CHARACTERISTICS OF DRAINAGE BASINS 


: By Bernard L. Golding, 1 M. ASCE and Dana E. Low,2 A. M. ASCE 


: SYNOPSIS 


the physical characteristics of major significance, such as; area of the 
mage-basin, slope of the principal drainage channel, length of the princi- 
fe inags channel, shape of the drainage-basin (basin shape), and general 

of the land in the basin (land slope), are reviewed and defined in terms 
a usage. A standard procedure for plotting the area-distance curves 
area-elevation curves is suggested. Also discussed is a new method de- 
ed for computation of a modified area-elevation curve. Numerical pro- 
res for computing the various important physical characteristics utilizing 
tandard area-distance curves and the modified area-elevation curves are 


uded. 


: INTRODUCTION 
enerally, the hydraulic engineer is lucky if there is astream-gaging station 
.e immediate vicinity upstream or downstream of the point in which he is 
ested. Thus in most cases, he is required to estimate stream flow by com- 

the drainage-basin above the point of interest with another basin having 
er physical characteristics. In the case of a highway, the hydraulic engi- 
‘ generally computes discharge corresponding to a particular frequency by 
ga graph for that frequency and entering with an argument in some way 


ste. Discussion open until August 1, 1960. To extend the closing date one month, 
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expressive of the physical characteristics of the drainage-basin. These gré 
have been developed by correlating discharge against physical characteris 
of many similar drainage-basins. In the case of a dam, the unit-hydrog! 
method is used. Coordinates of the peak value of the unit hydrographare tz 
froma graph in which the peak (maximum y-coordinate) and lag (correspon 
x-coordinate) of unit hydrographs for basins having similar physical cha 
teristics have beencorrelated with the physical characteristics of those bas 

Many recent papers have been published using such correlations so th 
knowledge of basin characteristics has become an essential for any engil 
working in the hydraulics field. 

These various physical characteristics of drainage-basins have been def 
in many different publications inthe past bymany different authors.3,4,9,6 ' 
has resulted insome of the same named basin characteristics being define 
several different ways. Also, many of these original publications are no k 
er obtainable as they are out of print. This has made the frequent referel 
to such publications in many recent papers quite useless. 

In an attempt to remedy this situation, the various drainage-basin cha: 
teristics are reviewed, and a suggested procedure for the computation of 
more important characteristics is presented. 


AREA OF THE DRAINAGE BASIN 


The area of the drainage-basin is its most important physical charac 
istic. The area is actually the horizontal projection of the land surface f 
which run off into the surface channels, above the point of interest, occurs 
is ordinarily measured on topographic maps by location of the divide, or s 
rating ridge, which determines whether surface run off will flow towar¢ 
basin in question, or flow into adjacent basins. The size of the drainage 
is generally expressed in square miles (or acres for small areas). Occas 
ally, sub-surface drainage-basins (ground-water flow) encompass diffe 
areas. However, this complication is very rare and generally not of ¢g 
significance, 


LENGTH OF THE PRINCIPAL DRAINAGE CHANNEL 


The length of the principal drainage channel is normally measured < 
the main channel, from the point of interest to the drainage-basin divide. 
length is normally measuredas short chords on United States Geological 
vey (USGS) quadrangle maps, which along with the projection of the prin 
channel to the basin divide, may cause some minor discrepancies in the l 
when measured by different individuals. However, such discrepancies 
generally insignificant in their effect on hydrologic computations. 


3 “Topographic Characteristics of Drainage Basins,” by W. B. Langbein, U. S. 
Survey Water Supply Paper 968-C, 1947. ; . 

4 “Drainage Basin Characteristics,” by R.E. Horton, Transactions Amer. Geox 
cal Union, Vol. 13, 1932, pp. 350-361. 

5 “Synthetic Unit Graphs,” by F, F. Snyder, Transactions Amer. Geophysical 1 
Vol. 19, Part 1, 1938, pp. 447-454, 

“Unit-Hydrograph Lag and Peak Flow Related to Basin Characteristics,” by 

Taylor and H. E. Schwarz, Transactions Amer. Geophysical Union, Vol. 33, 195 
235-246. : 
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SLOPE OF THE PRINCIPAL DRAINAGE CHANNEL 


The slope of the principal drainage channel is probably the most significant 
7sical characteristic after the area of the drainage-basin. The simplest 
thod for expressing the slope of the principal drainage channel is to divide 
length of the channel by the difference in elevation between its upper and 
yer ends (Definition 1, Fig. 1). Generally, the length of the channel is mea- 
-ed beyond the upper end of the clearly discernible stream channel to the 
uinage divide and then divided by the difference in elevation between this 
nt on the ridge line and the point of interest on the channel. The USGS has 
2d a slope parameter instudies inthe New England states in which the length 
that part of the stream between points 85% and 10% of the total distance above 
point of interest is divided by the difference in elevation at these places.” 
ey are of the opinion that the most upstream part of the slope, in the steep 


ELEVATION - DISTANCE CURVEX. 


SLOPE BY DEFINITION I~ 
SLOPE BY DEFINITION 2 
SLOPE BY DEFINITION 3 
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FIG. 1.—SLOPE OF THE PRINCIPAL DRAINAGE CHANNEL 


adwaters, might affect the slope out of proportion to the volume of water fur- 
shed by the headwater area, and also that the most downstream part of the 
ype of the principal channel may be so flat that it might not indicate the true 
rtion of the flow which the lower portion of the drainage area actually does 
7. ee he the principal drainage channel canalso be determined by dividing 
ice the mean height of the channel profile above the point of interest by the 
incipal channel. ; 

| aerenanens of the profile is obtained by integration of the Seete 
stance curve (channel profile) and dividing by the length pt the channel le if 

‘ion 3, Fig. 1). The integration can be accomplished by planimetering” e 
ea below the elevation-distance curve. This slope is also given by the slope 


| i ! M. S. Benson, Proceed- 
7 “Channel-Slope Factor in Flood Frequency Analyses,” by M. 58. a. 

s Paper 1994, pournal of the Hydr. Div., ASCE, Vol. 85, No. HY 4, April, 1959, pp. 
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of the line drawn through the origin of the stream-profile curve such that tt 
area under it equals the area under the profile curve. The mean channel heigl 
is at a pointon the above-cited line at one-half the length of the principal drair 
age channel measured from either end. Rather than computing mean chann 
height as such for use in determining channel slope, the mean basin height ( 
be defined, subsequently, under the heading “Land Slope”) is often used and tt 
resultant slope is called the mean basin slope. This effectively introduces 
weighting factor to account for slope of the basin aswell as slope of the strea: 
channel. 

A probable better and more rational definition for slope of the princip 
drainage channel is the slope of an equivalent stream having the same trav: 
time and same length (Definition 2, Fig. 1). This slope is computed by 

Sst = et Ph ne ah cs ee 
1 


Zp si 


where p equals the number of equal reaches into which the stream has bet 
divided (often 10), sj is the average slope of each such reach measured as tl 


change of elevation over the reach divided by its length, and Xp - represen 
i 
the summation of the Ve -values for each reach. Eq. 1 is based on the as 
i 

sumption that velocity is proportional to the square root of the slope. Th 
method of determining slope was defined by A. B. Taylor and H. E. Schwar 
M. ASCE and is used by the Corps of Engineers.8 

The latter two methods of slope determination give values thatare approx 
mately equal and which may be used interchangeably for small rivers (Fig. 1 


BASIN SHAPE 


The shape of a drainage-basin can be expressed in terms of various eas! 
defined and measured factors. One measure of basin shape is given by t 
quotient of the drainage area divided by the length of the major watercours 
A/L. Another measure of basin shape is given by the length of the major wate 
course squared divided by the basin area L2/A. 

The most commonly used basin shape factors are Leg and the product 
times Leg to some power, usually (L Leg)9-30, where L is the length along t 
longest watercourse and Lea is often defined as the distance along the m: 
drainage channel from the point of interest to a point opposite the comput 
center of gravity (centroid) of the drainage area.9,6 A more accurate definiti 
of Leg is given subsequently. The term Leg is computed by integration of 1 
area-distance curve and dividing by the drainage area. 

To construct the area-distance curve (Fig. 2) the channel lengthis first st 
divided into separate reaches. These reaches may be of equal length along 1 
watercourse, such as 1-mile reaches or equal parts of the total length. Div: 
ing points along the stream couldalso be made at points where contours cro: 
A preferable method, especially for larger rivers, is to divide the watercou1 
length just above and below each major tributary. The areas in this case < 
the contributary drainage areas at these points. After the channel is divic 
into reaches, the sub-drainage basins contributing to the stream between | 
limits of each reach are determined. | 


8 Unit Hydrograph Compilations, Project CW 153, Washington District, Corp: 
Engrs., Washington, D. C., 3 Volumes 1949, 1 Volume 1954. 
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A curve of channel distance upstream from the point of interest (or per- 
tage of channel distance) versus area (or percentage of area) contributing 
hin that distance is then drawn as shown in Fig. 2. The distance to the cen- 
id of the area is found by dividing the area above the curve by the drainage 
a (or percentage of the drainage area) which is the maximum ordinate. As 
alternative, Lea may also be found by integrating the area-distance curve 
dividing by the drainage area. The integration is accomplished by “plani- 
tering” the area above the area-distance curve. 

A suggested method for plotting the area-distance curve and a sample com- 
ation of L,, based on this curve will be presented subsequently. 

Physical Definition of Lcg-—Previously, it was stated that the centroid of 
area above the area-distance curve determines Log. In other words, Leg 
the length measured up the stream channel from the base of the drainage 
ato a point corresponding to the centroid of the area above the area-distance 
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DISTANCE UP PRINCIPAL CHANNEL FROM BASE 
FIG. 2.—BASIN SHAPE FACTOR, Loa 


ve. Although all of the various methods of approach that utilize the Lea 
stion derive it in the same general manner, some writers have added physi- 
definitions of Leg which are not rigidly correct while others leave the 
sical interpretation of Leg purposely vague. For example, Lea has been 
ined as “the distance along the stream channel to a point opposite the cen- 
idof the drainage area.”6 This definition is actually incorrect since the lo- 
ion of the centroid of an area depends solely on its shape, and Lea, derived 
m the area-distance curve, obviously depends on other factors as well. oe 
other hand, L,.g is sometimes vaguely definedas “the distance to the center 
: ” 
Physically speaking, just whatis Leg? Fig. 3 shows a hypothetical Ree 
a. Axis x-x is drawn perpendicular to the major axis of the area throug 


) « Applied Hydrology,” by R. K. Linsley, M. A. Kohler, and J. L. H- Paulhus, 
3raw-Hill, New York, 1949, p. 456. 


oe - 


6 March, 1960 HY 


the base of the area. The true area-centroid is then located a distance ¥ abov 
this axis where 


in which AA is an elemental area, yj is its distance from the axis x-x, and. 
is the total area. The term Leg can be expressed by a similar formula 


in Eq. 3 yg is the distance measured along the stream to the point on the chan 
nel at which overland runoff resulting from rainfall on AA first reaches th 
stream. 


! 
| -strean Channel 
a 


Direction of runoff from A A 


FIG, 3.—HYPOTHETICAL DRAINAGE AREA 


LAND SLOPE 


Defining a characteristic which adequately describes land slope has prove 
difficult. However, this characteristic is a major factor in indicating the trave 
time of overland flow, so several methods for determining land slope have bee 
developed. 

Most of the methods generally used for defining land slope take into accoul 
the general slope of the entire basin including the tributary stream netwo1 
contributing to the major watercourse, but excluding the major watercourse, 

Average land slope has been defined as equal to the plotted contour interv: 
times the total length along all contours in the drainage-basin divided by tt 
drainage area. However, measuring the length of contours proved a tediot 


and time-consuming operation and hence the following short-cut methods wel 
derived and are used: 
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In the intersection-line method, a grid of uniformly spaced parallel and 
‘pendicular lines is laid over a contour map of the drainage area. The num- 
>of contours crossing each subdivisionof the grid withinthe areais counted. 
e land slope in either grid direction is then s = uae 
nber of contour crossings forall lines in one direction, AZ denotes the con- 
r interval, and L is the total length of grid lines in one direction within the 
2a. The average land slope is then given by some relation between these two 
ues of slope or by 


where N is the total 


N' AZ 
Li 


s'= 


2re N' and L' are the sum in both grid directions of the values previously 
inedas N and L, and @ is the averageangle of intersection between the con- 
rs and grid lines. Often a value of 1.57 is used for sec @ as it is the av- 
ge secant of angles from 0° to 90°. The Corps of Engineersuses this method 
‘inding land slope.8 


ELEVATION 


DRAINAGE AREA ABOVE 


Mean Elevation 
(Mean Basin Height) ] 


FIG. 4.—LAND SLOPE 


Another method of defining land slope, or more accurately, basin slope, util- 
3 an area-elevation curve in which the drainage area above a certain ele- 
on is plotted against that elevation. Area-elevation data can be assembled 
planimetering the area enclosed by each contour and the basin divide and 
ted as shown on Fig. 4. 

[he mean basin height can be obtained by integration of the area-elevation 
ve and dividing the result by the drainage area. It can also be obtained by 
method shown in Fig. 4 in which a line is drawn such that AY = Ag + Ag. 
intersection of this line with a line drawn parallel to the abscissa at 50% of 
drainage area indicates the mean basin height. Mean basin slope is then 
nedas two times the mean basin height divided by the length of the channel, 
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Another curve can be plotted conveniently by using the same data derive 
for producing an area-distance curve. This is called a modified area-elevatic 
curve. Rather thancutting the drainage-basin into slices between adjacent cor 
tours, it can be divided into sub-drainage basins contributing to equal reache 
along the length of the principal drainage channel to correspond with the methe 
used in computing the area-distance curve. The preferable method, especial) 
for larger rivers, is to break the channelat the tributaries just above the poi 
where they jointhe main channelas described under the heading “ Basin Shape. 
In this case, area no longer refers to area lying above a certain elevation b 
rather to the area contributing runoff to the major stream above the point : 
which the stream bed is at a certain elevation. Thus the mean basin height - 
obtained actually weighted by distribution of drainage area long the stream prec 
file. Mean elevation and mean basin slope are correspondingly weighted. 

The foregoing method saves considerable computation time over the usu 
method for computing an area-elevation curve by utilizing one set of area da 
as the basis for both curves. In the opinion of the writers this results in 
more meaningful mean basin height and mean basin slope. It is undoubted 
not applicable in studies of snow hydrology where the temperature-elevatic 
curve makes true land elevations more significant. Further study may als 
indicate that the proposed modification should not be used in certain other sul 
divisions of the subject of hydrology. 


SUGGESTED PROCEDURE FOR DEVELOPING CURVES 


In the plotting of the area-elevation curve to compute mean elevation ai 
mean basin height and the area-distance curve to compute leq, it makes | 
difference whether zero channel distance is taken at the base of the draina; 
area or at the drainage divide. Also it does not really matter whether t 
drainage area above a certain point on the channel or the drainage area belc 
that point is plotted versus that channel distance. The plotting units may | 
actual or percentage values and the location and sense of scales could co 
ceivably vary in any one of several ways. However, in order to eliminate ai 
possible confusion a suggested standard procedure similar to the method us 
by the Corps of Engineers8 is outlined herewith: 

Channel distance will always be measured up the channel fromthe base a 
will be expressed in percentage of total length. Elevation will always be me. 
sured above the stream-bed elevation at the base of the drainage-basin and w 
be expressed as the percentage of the total elevation drop from one end of t 
stream to the other. Drainage area will also be expressed as a percentage 
the total. It will be recalled that in the development of the area-distance r 
lationship the term area referred tois the area of the drainage-basin contribu 
ing below the point along the stream in question. On the other hand the are 
elevation relationship required the consideration of the contributary area abo 
the point in question. This distinction need cause no confusion in actual pra 
tice, as will be shown. 

The fact thatall scales go from 0% to 100% enables the engineer to plot be 
curves on a square piece of graph paper. Since both curves appear roughly 
diagonals on such a graph, the scales are oriented so that they will cross or 
once, thus eliminating any danger of confusion (Fig. 5). The value of lea 
thenfound by measuring the area above the area-distance curve in square in 
es and converting to square percentage, for example, for horizontal and ver 
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| Scales of 1 in. = 10%, the conversion factor is 1 sq in. = (100%)2. This 
ure divided by 100% is the Leg expressedas a percentage of the total length. 
To find mean basin height (mean elevation), the area above the area-elevation 
tye is determined, converted to square percentage and divided by 100%. The 
sult is the mean basin height expressed as a percentage of the total height. 
‘an basin slope is then mean basin height divided by one-half the total chan- 
_ length, both values having been converted from percentage to actual lengths. 
The modified area-elevation curve is plottedin a similar fashion to that de- 
sibed for the standard area-elevation relationship. 


% ELEVATION 
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MODIFIED AREA-ELEVATION CURVE 
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FIG, 5.-HYPOTHETICAL DRAINAGE BASIN 


a modified area-elevation curve 


E _ ints for plotting 
ees aa tor. 1 drainage-basin are shown in 


-distance curve for a hypothetica 
3 1 " Ae on the stream were made ateven tenth aha along oa 
ilchannel length. Other methods, as previously noted, could have anes 
e that distance is measured in anupstream direction and that are ep Lip 
ion are measured below the point on the channel, The ‘4 nes = 
vation curve is plotted using the lower right-hand corner as the origin. 
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An alternate solution is to add another column to Table 1 for percentag 
area above, which is determined by subtracting percentage area below fron 
100%. The modified area-elevation curve can then be plotted directly using th 
origin in the upper-right corner. 


TABLE 1.—HYPOTHETICAL DRAINAGE BASIN (DRAINAGE AREA = 8.10 SQ MILES, 
CHANNEL LENGTH TO BASIN DIVIDE = 5.42 MILES) 


Length, L, % L Elevation, E, %E Area, A, in %A 
in feet in feet below square miles below 
(1) (2) (3) (4) (5) (6) 
0 0 21 0 0 0 
2,860 10 43 11 0.70 8.6 
5,720 20 54 16.5 2.43 30.0 
8,590 30 63 21 3.14 38.8 
11,450 40 78 28.5 5.07 62.6 
14,310 50 94 36.5 5.28 65.2 
17,170 60 116 47.5 5.62 69.4 
20,030 70 134 56.5 5.84 Took 
22,900 80 180 79.5 6.55 80.8 
25,760 90 202 90.5 T.51 O227. 
28,620 100 221 100 8.10 100.0 
TABLE 2 


Elevation, 


Reach length, 


in feet 


(3) 


in feet 


0.00770 0.0878 


54 0.00385 0.0621 16.1 
63 0.00315 0.0561 17.8 
78 0.00525 0.0725 13.8 
94 0.00560 0.0748 13.4 
116 0.00770 0.0878 11.4 
134 0.00630 1 
180 0.01608 1G 
0.00770 1. 


0.00735 


A sample computation of the weighted values of mean basin height, me: 


elevation and land slope (mean basin slope) based on the data given in Table 
is as follows: ; 
Area above modified area-elevation curve = (3810%)2 
Mean Basin Height = (0.381) x (221 - 21) = '76.2 ft =e 
Mean Elevation = 21 + '76.2 = 97.2 ft (m.s.1.) 
Land Slope (Mean Basin Slope) = 76.2/14,310 = 0.00532 ft per ft 

= 28.1 ft per mile 


DRAINAGE BASINS 11 
1¢ computation of Leg is as follows: 


Area above area-distance curve = (4330%)2 
Leg = (0.433) x (5.42) = 2.34 miles 


atermination of the Channel Slope.—In Table 2 stream slope for the same 
thetical drainage basin considered in Table 1 is computed by the method 
iously reviewed in which the slope of the principal drainage channel is ex- 
sed as the slope of an equivalent stream having the same travel time and 
length.6,8 Pertinent information from Table 1 is thus repeatedin Table 2. 


SUMMARY 


review of the more important physical characteristics of drainage-basins 
een presented along with a suggested standard procedure of plotting the 
rtant area-distance curve, from which Leg canbe computed, and the area- 
tion curve from which the land slope (mean basin slope) can be computed. 
ditional characteristic, developed fromthe modified area-elevation curve, 
een described. 
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DEVELOPMENT OF FLOW IN TANK DRAINING 


By David Burgreen! 


SYNOPSIS 


theoretical study is made of the development of velocity with time when 
ve at the end of a drain pipe is suddenly opened. The effect of friction, 
surization, and length of drain pipe is examined. The true velocity, ob- 
d by considering the fluid inertia, is compared to the velocity obtained by 
nore well-known method which disregards the fluid inertia and assumes 
the draining is governed only by the prevailing head. 


INTRODUCTION 


1e problem treated herein is a generalization of a practical problem 
ng with the shutdown of a reactor by the rapid draining of moderator or 
ctor liquid. It is important to know quantitatively the rate of draining of 
iquid from the tank to ascertain that it is sufficiently high to give rapid 
or shutdown. The mechanics of rapid liquid-level change or of rapid tank 
ing falls in the realm of non-steady flow. The moderator or reflector 
lis initially at rest, and the order of magnitude of time required to reach 
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peak velocity is the same order of magnitude of time, from a nuclear point 
view, in which the shutdown must be completed. 

The more well-known type of tank-draining problem is the one in which 1 
rate of discharge is governed by the decreasing head inthe tank. Althot 
the rate of flow varies with time, this is not a true non-steady flow proble 
but rather one of slowly varying steady flow. In flow under decreasing he 
the inertia forces are neglected, whereas in non-steady flow the presence 
inertia forces is implied. The flow relationships that are derived without 
inclusion of inertia forces are designated as “terminal flow” relationships 
this paper, since they represent an approximation of the terminal phase of 
flow. 


H = initial height of 
liquid in tank 


x | 


L, = vertical length 
of drain pipe 


Pog 
i} \ 
3 
Horizontal length 
of drain pipe 


FIG. 1.—TANK AND DRAIN PIPE ARRANGEMENT 


A related problem dealing with the hydrodynamics of a conservative n 
and energy system has already been studied.2 The present problem 
with a non-conservative mass and energy system, in which mass and en 
in the form of draining liquid is continuously being lost from the system. 


BASIC EQUATIONS 


The system consists of a large pressurized tank from which liquid dz 
through a drain pipe of smaller diameter having a vertical and a horiz 
length. It drains under the influence of gravity, assisted in some cases 
gas pressure on the surface of the liquid. Fig. 1 shows the arrangeme 


2 “Hydrodynamics of Liquid Poison Scram System,” by David Burgreen, Nu 
Science and Engineering, Vol. 4, No. 1, July 1958. | 
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ind drain pipe. The summation of forces on the slug of liquid in the tank 


Ww d2y 
1 ie Py Ay+w (H-y)A, a (H-y) Ai a2 MONS Sieivice oa nus cones, oot 6 (1) 


is assumed in deriving Eq. 1 that friction forces retarding the flow in 
ink may be neglected. A balance of forces on the slugs of liquid in the 
sal and horizontal lengths of drain pipe give, respectively: 


Ww d2x ow Ly /dx 7 

eee F 2725 Sa lyn hy Ag so 28 la) =0.. (2a) 
w d@x ow Lg /dx\2 

Meme 0279. 51242 9-2 ¢ sot 5 () OR ete) <0), e 8) 61 eke (2b) 


ich x represents the displacement of a fluid particle in the drain pipe 
3 most clearly defined by the continuity relationship A, y = Agx. The 
on of Eqs. 2 yields 


, w. dx  w dx ¥ 
a Se Ss Cee rc} so (2) = (0) 5 ceteris era Goatees 45 (3) 
j- 3, L=1Ly + Lo and F represents the number of velocity heads lost 
th pipe friction and would include entrance, exit, bend, and valve losses 


lition to the straight pipe losses. At the junction of the pipe and tank the 
ure-velocity relationship is the same as for steady flow.2 It is 


= Py “35 (=) [1- (2 ee AE Ee 5 « tae lett (4) 


ded the length of the drain pipe is greater than its diameter. 
e use of the continuity relationship, Ay y = Ag x, together with Eqs. 1, 3, 
produces: 


Ag (A2\2_ Tax 1 Ag\2] /ax \2 
mone | )>x|—~ +5 | 4F-(-- ss 
Ay Ay at2 2 Ay dt 
Pp 
-¢ [ae Phe ty - 2 x | PA Vir by. ee eceCr Ee ih Ce ONO (5) 
Ww Ay 


. 5 is expressed in terms of the displacement of a particle in the drain 
The flow could be expressed just as readily in terms of the drop of level 
tank, y- In deriving Eq. 5, Poo, the pressure at the end of the drain pipe 


een assumed to be atmospheric (zero gage). 


INITIAL ACCELERATION 


the start of flow the displacement and the velocity are zero. From 
the initial acceleration is 


Bee eae Peeper aati 887 ee * 8 (6) 
at? / ¢-0 L+RH 
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where R = Ag/A, and n H = Po1/w. The pressurization in the tank is thus 
pressed as the number, n, of initial tank heads. The initial acceleration of 
liquid in the tank, which is the same as the acceleration of the dropping s 
face of liquid is 


2 
(<3) Are boar ed Pes een 
dt2 /;_9 L+RH 


It is clear from Eq. 7 that without pressurization the initial acceleratior 
the liquid in a tank, to which a drain pipe is attached, can never be as higk 
g, since R <1. In order to obtain an initial acceleration of g, the requi 
pressurization n H is 


*) Lg + (1-R) Lj 
= R 


a a ee 8 eke a 0.0 Bye CR 6. oe 1 CC _S eee 


CROSS-SECTIONAL AREA OF A VERY SMALL DRAIN PIPE 


When the flow area of the drain pipe is very small in comparison to 
flow area of the tank, Eq. 5 simplifies to 
d2x  14+F ‘Gah | 
ee te dt - g[ M+1) H+ Ly] =0 ocee ue ss oe 
and represents flow under constant head with inertia effects. It is thus indi 
tive of the initial flow behavior only, since the approximation of Eq. 9 is 
valid for large x. The first integrals of Eq. 9 with respect to time and¢ 
placement are, respectively, 


1/2 
. poe a+ tal) tanh [ {2 ef (n+l) H+Ly] (1+? unt 2 


dt 1+F 2L 

and 

dx [2 e[ (+1) H+Ly] 71/2 -[ (1+F)/L 1/2 

oe pe See aka eee Beech Ate ee 


and the displacement may be written explicitly as a function of time by ct 
bining Eqs. 10: 


-1 
xe i In | 1 - tanh2 ({2 g[ (n+1)H+L ] (1+F)}1/2 ‘<) | = pica Oa a 


d+ 


If there were a horizontal drain pipe only and no pressurization or ee 
then Eq. 10a would become 


Se = [2¢ Hea]? tanh | @ gh? St a ee 

2 Lo Lo ; 

Eq. 12 has been one from the Bernoulli equation for non-steady flow. 

From Eq. 10b it is seen that for small values of x/L the velocity 7 
drain pipe may be approximated by 


ax zl 2 g x[ (n+1)H+L] | 1/2 
dts ve eae eee 


ee a a eh LR ee a 
3 «Fundamentals of Hydro and Aeromechanics,” by Prandtl-Tietjens, McGraw 1] 
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result is compatible with Eqs. 6 and 7 and is the velocity associated with 
stant initial acceleration. When the drain pipe area is extremely small, 
hange in head is very slow and the starting velocity of the terminal flow 
t sensitive to fluid movement in the lower drain pipe, and can be ex- 


ed as 
ee ee ram sse eso 4 


TERMINAL FLOW 


e solution of Eq. 5 with the inertia term d2x/at2 = 0 is the solution of the 
known tank-draining problem in which the discharge is governed by the 
sntary head and the inertia of the fluid to a change in velocity is disre- 
d. It is of considerable practical importance and, as will be shown, it is 
y close approximation of the flow for the case when both the area ratio 
all and the length of drain pipe is small. Under zero inertia conditions, 
lution of Eq. 5 in terms of the y-coordinate, or liquid level in the tank, is 

dy _ E g[ (n+1)H+L4-y] | 1/2 

dt 1+F-R2 


ntegral of Eq. 15 gives the time required for a specified drop in liquid 
It is 


1/2 
|? (HF-RY) | é | (a+) 4-¥0) 1/2 5 (ns1)H+Ly-y]!/2 |... (16) 


> Yo is the liquid level at t = 0. 


GENERAL SOLUTION 


though Eq. 5 is a non-linear differential equation with variable coeffi- 
3, the first integral, in terms of the tank level y, is readily obtained. It 
the velocity as 


2 
-| ZeR° [ msayeenys OED | [1- (2- ah 
14+F-R2 14+F-2 R2 


) ores 


PreRy sie Seah UUER,, eae Sek ALES Ber PCTS (17) 
14+F-2 R2 


7 has four parameters in addition to the independent variable y/. H. The 
of these parameters will be examined separately. Consider first the 
on when the drain pipe lengths Ly = Lg = H, when there is no pressurl- 
1, and the friction forces are neglected. Eq. 17 becomes 


2 
‘R2(24+2 R-3 R2 y R \(1/R%)-1]_y x eet 
=(1-R2)(1-2 R2) [2- (1- H 2+R ) | H 1-2 R? an) 


i i in pi to tank area. 
3 is plotted in Fig. 2 for several ratios of drain pipe area 
this ratio is unity, the configuration represents a pipe elbow of uniform 


4 
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diameter whose vertical length is twice the horizontal length. Eq. 18 t 
more simply expressed as 


v2 y y 
eS ee - —— A: Pe ee 
2gH E+in(1 ne 


A more general expression derived from Eq. 17, neglecting friction and pr 
surization, for the discharge of fluid from a pipe elbow of uniform diam« 
of total length L + H and horizontal length Lg = L - Ly is 


ova oy. L y H 

H L+H /’ 
In Fig. 2, a free-fall curve, which is included for comparative purpos 
shows the velocity of a free-falling body, dropped from a height H, at any € 


Rate of Drop of Liquid Level 


on Relative Liquid Level 


FIG, 2.—VARIATION OF FREE SURFACE VELOCITY WITH LIQUID 
LEVEL IN TANK. Ly =Lo=H: n=); F =0. 


vation. It is a special case of Eq. 20, or the R = 1 curve when Lg, the h 
zontal length of pipe, is equal to zero. The curves in Fig. 2 show the rat 
fall of the liquid surface in the tank to decrease with decreasing area rat 
The range in which the flow is governed primarily by inertia forces also. 
comes smaller with decreasing area ratio. This may be observed meee 
increasingly early blending of the true flow curves with the terminal 
curves as the area ratio becomes smaller. The dotted terminal flow cul 
are plotted from Eq. 15 which neglects inertia effects. These are fair 
proximations of the terminal stages of the flow. 
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en the area ratio is quite large as, for example, when R = 1/2 the 
ity-displacement curve shows that the inertia forces primarily deter- 
the character of the flow. For large area ratios the velocity- 
cement curves resemble the free-fall curve and any approximation with 
nal flow would not be valid. Note that the true velocity-displacement 
s indicate, properly, that the velocity starts at zero and builds up, while 
srminal flow curves indicate, incorrectly, that at the start of flow the 


ity is V = (2 g H)1/2. 


= Rate of Drop of Liquid Level 


_vy 

¥ 2¢H 
u 

las 


* 
| 


+ = Relative Liquid Level 


_ FIG. 3.—VARIATION OF FREE SURFACE VELOCITY WITH LIQUID 
LEVEL IN TANK; PRESSURIZATION: n=10; F=0; Ly =Lg=H. 


e effect of pressurization is shown in Fig. 3. The pressure on surface 
liquid is assumed to be constant and is arbitrarily taken as ten tank 
that is, n= 10. These curves show that a rate of drop in liquid level, 
ess of the free-fall rate is easily obtained. It will be observed that for 
lected pressurization, an area ratio of 1/4 produces a rate of drop in 


_ 
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liquid level approximating free fall. The terminal flow approximation is | 
rather poor even for small area ratios. 

The result of taking friction into account is demonstrated in Fig. 4. 
total friction loss in the system is arbitrarily taken at three velocity he 
(drain-pipe velocity). Comparing the curves of Fig. 4 with those of Fig: 
and 3 it will be noted that the velocities, for corresponding area ratios, 
smaller and that terminal flow is established sooner. Although terminal 1 
is now a better approximation to true flow, it is obviously not a satisfact 
approximation, for moderate area ratios, of the early stage flow. 

The effect of drain pipe length on the development of flow is well dem 
strated by considering a configuration having a small area ratio typical | 


v 
Tia Rate of Drop of Liquid Leve 


Oo 6860 0.2) 603 04 OLS «(0.6 0.76 (0:85 0.9 =1.0 


* = Relative Liquid Level 


FIG, 4.-—VARIATION OF FREE SURFACE VELOCITY WITH 
LIQUID LEVEL IN TANK, Lj= Lo =H; n=0; F=3, 


water tank with drain pipe attachment. An area ratio of 0.01 is selected, 
the velocity-displacement relationship is plotted in Fig. 5 for lengths of d 
pipe of zero, H, and 5H. The curves show, in common, that the ine 
transient is completed and terminal flow started when the maximum rat 
flow, more or less, has been attained. } 

The curve in Fig. 5a represents the development of flow when a very. 
drain pipe, equal to 5 H, is attached to the tank. The peak velocity is not 
tained until about half the tank has been drained. This indicates that a 
drain pipe produces a long inertia transient, and that computations of i 
charge based on the prevailing head lead to incorrect results — even wher 
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ratio is small. When the drain pipe is of moderate length, equal to H, 
the peak velocity is reached when one-tenth the tank has been drained. 
is shown in Fig. 5b. 

@ fundamental differences between inertia flow and terminal flow are 
ht out by the curves in Fig. 6 in which the fluid acceleration is plotted 
st the drop in liquid level. Since terminal flow is flow without accelera- 
the deviation from terminal flow is measured by the magnitude of the 
eration. The parameters are the same as those used in the velocity- 
ywcement curves of Fig. 2; namely, unpressurized frictionless flow with 
Lg = H. The uppermost curve is a straight horizontal line representing 


a 


TH 


x = Relative Tank Level 


FIG. 5.-DEVELOPMENT OF FREE SURFACE VELOCITY WITH 
LIQUID LEVEL IN TANK. R= 1/100;n=0; F= 0. 


al the equation (1/g)(d2 /dt2) = 1. An examination of the curves will 
Biscsy a case ue Sea 9 is a maximum at the start of flow 
hereafter decreases continuously with increasing displacement. As the 
ratio becomes small the curves show a range of decelerating flow which 
sponds to the region of decreasing velocity. Terminal flow a ee 
d by the abscissa, (1/g) (d2y/dt2) = 0. The general trend of the eu vgs 
yard zero acceleration as the area ratio approaches zero. This ees 
_ will prevail also for configurations other than Ly = L2 = H, but w. 


r drain pipes will tend toward zero acceleration at a slower rate. 


s 
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To obtain time-displacement curves it is necessary to integrate Eq- 
Analytically, it appears to be possible to perform the integration only for s 
cific values of F and R. In general the time- displacement curve would be 
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FIG, 6.—VARIATION OF ACCELERATION OF FREE SURFACE WITH 
LIQUID LEVEL IN TANK: Lj =L =H; n=0; F=0. 


tained by numerical integration of the velocity-displacement curve. The t 


the reach a given level, yp, is 


Initially 1/V tends to infinity. This represents a difficulty for numerical 
tegration. For small values of y/H Eq. 17 becomes 


y2 ~2ey Ri (n+1)H+L1] 
ae 


L+R H 


An examination of Eqs. 6 and 7 indicates that Eq. 22 is simply 


sats 


H 
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he expression for velocity when the acceleration is constant. This ap- 
imation may be used to start the integration, since for constant accelera- 


2Qy /2 


dy 
dt2 
small values of t 


mn 


DISCHARGE THROUGH A WEIR ANNULUS 


a 
t is fairly apparent from the foregoing material that to obtain the rapid 
nage required for the rapid shutting down of a reactor it is necessary to 


| | y = surface level drop 


Weir drain annulus 


Liquid surface 


FIG. 7.—CONFIGURATION OF WEIR DRAIN SYSTEM 


a large discharge area. Rather than use a very large pipe or group of 

3, it is expedient to use an annular weir discharge as shown in Fig. 7. In 
arrangement the length of the drain pipe L = Ly + Lg is assumed very 
1 in comparison to (Ag/Aj)H. The solution of Eq. 5, disregarding pres- 
zation and friction, is 


é 2)-1 
Me R? y J ras | Lay SNE REE 5 25 
aie wise Re [i-#- (2) ke ie) 


45 is plotted in Fig. 8 for ratios of drain pipe area to tank area of re 
1/8. The corresponding terminal-flow curves are also shown. Note tha 


’ 
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now, in the absence of any length of drain pipe, terminal flow is a better aj 
proximation of the true flow, and for a practical application, would give 
good approximation of the time required to completely drain a tank. It i 
however, not satisfactory when the interest is primarily in the starting flo 
which it does not approximate very well. 

Eq. 25 can be integrated to give the drain time explicitly. When R < 1/1 


the drain time is 
2 
R 1 
(saat) 
eG 6 Mie V2 1 4 Re) a 
2g 1-2 R2 ee ) 
2 
in which Jis the gamma function. WhenR > 1 /V2 


5 ( a R2 ee (22 )r() pe (2¢ 


2 R2-1 2-5 R2 
chen 


For R = 1/V2, Eq. 25 becomes 


v2 
MENS = - ye oe 0 
= (2 #) In oat Wee ties greet Pe ( 


As 1/V¥2 is a very large area ratio, it is well represented in the early ste 
flow by free-fall. This is shown in Fig. 9. 
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FIG. 8.—-VARIATION OF FREE SURFACE VELOCITY WITH LIQUID 
LEVEL IN TANK, L<<(A2/A,)H; n=0; F=0, 


3 
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The time required for the complete draining of the tank when R = 1/V2 
obtained by integrating Eq. 27 or from Eqs. 26. It is 
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10 shows the elevation of the liquid surface at any time. It is plotted for 
1/V2 , for free-fall, and for (improperly assumed) terminal flow. Note 
until about three tenths of the tank has been drained the time-displacement 
tionship is practically that of free-fall. The divergence becomes more 
ificant as the draining of the tank nears completion. To complete the 
ning, 25% more time is required than the time it takes for a body to fall 
ly through this distance. Terminal flow is not at all representative of the 
state of affairs. The fact that the free-fall and terminal flow curve meet 
/H = 1 is fortuitous. 


Rate of Drop of Liquid Level 


Be eN 


0. 
Relative Liquid Level 


FIG. 9.—VARIATION OF FREE SURFACE VELOCITY Re as 
LEVEL IN TANK. L<<(A2/A,)H; n=0; F=0; R=1/V2. 


CAVITATION 


t the uppermost point of the drain pipe, near the junction with the ae 
pressure in the liquid may be quite low at some stage of eae - 
sure at this point is less than the vapor pressure of the liqui , =e as 
ration or cavitation may occur. The pressure at this point is 0 aine 


1 Eq. 3. It is 


a \2 
P2_ Po2 | _b dx + (F) owe 6 <a Bo Bee. 6 '@ eKehe (29) 
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In Eq. 29, Pog is atmospheric pressure and Po is the absolute pressure. 
P., is the vapor pressure of the liquid then (Pg9/w) - (P.;,/w) = her wot 
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FIG, 10.—CHANGE IN LIQUID LEVEL WITH TIME, L<< (Ap/A,)H; 
n=0; F=0; R=1/y2. 


represent the height of a column of liquid that can be supported when separ: 


tion is about to occur. Separation will not occur as long as (Po,/w) < (Po/ 
or 


P Pp 2 2 
ee ee ee 2 ae ( 

Ww Ww 1g. geal v2e@ \V dt ; 

and separation will occur when 
y 

stig K cy Ht EN” 

1 g at2 2¢ dt Leh te o's 6" Se a 2 8 we a ee . 6 

a 


Thus when the fluid is stationary and not accelerating (d2x/dt? = dx/dt = 
separation occurs when the vertical leg Ly exceeds hey. At the start of flor 
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n dx/dt = 0, in order for separation to occur, it would be necessary that 


L d@x 
ee err ai Re eo Ses Shee kee (32) 
\t any given time or at any liquid-level displacement there will be an as- 
ated acceleration and velocity. When the acceleration and velocity fail to 
sfy the inequality, Eq. 31, it will demonstrate that no separation has taken 
e. When separation takes place the draining of the tank is no longer 
srned by Eq. 5. 

3efore the pressure at any point in the pipe reaches the absolute vapor 
ssure of the fluid and produces flow separation, it will pass through a pres- 
> range less than atmospheric. There would then be a tendency for air to 
< up through the pipes. However, as in the case of a siphon, as long as 
e is a substantial flow velocity in the drain pipe the air will not move 
iter the flow and the pipe will flow full. 


SUMMARY 


‘he major findings of the foregoing study of tank draining are the follow- 


. When the ratio of drain pipe area to tank cross-sectional area is small 
the length of drain pipe is short, then the conventional manner of comput- 
the rate of draining, which disregards the inertia of the fluid, may be used 
otain reasonably accurate results. 

. When the area of the drain pipe is a large fraction of the tank area, the 
of draining is given by Eq. 17. The neglect of the fluid inertia, in this 
», leads to incorrect results. 

. When the length of the drain pipe is a large fraction of, or greater than 
initial height of fluid in the tank, it is necessary to use Eq. 17 to compute 
rate of draining, and neglect of the inertia of the fluid, here again, leads 
correct results. 

. Since either large area ratios or long drain pipes have a strong influ- 
on the flow, a combination of moderate area ratio and drain pipe length 
requires Eq. 17 for the proper computation of the rate of draining. 

. As one may expect, the effect of drain-pipe friction is to decrease the 
of draining while the effect of tank pressurization is to increase the rate 
raining. For expediency the form of the friction term was taken as the 
dy state form, although the flow is actually non-steady. 


APPENDIX. —DERIVATION OF MOMENTUM EQUATION BY 
ENERGY CONSIDERATIONS 


fhe kinetic energy of the liquid in the tank is 


w Ay (H-y)¥? Nie 


K.E. = 2¢ 
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and in the drain pipe it is 


roe w Ag (L4+L2) x2 _wAgL x2 (3: 
ire 2e ermal oe eek , 
The potential energy, using x = 0 as a datum is 
P.E. = 5 Ay (Hey dies ait dlnige- 20 ot cerieenale oem (3: 
for the tank liquid and 
Ww 2 
P.E. = “3 Ag L, - WAg Ly lg -- eee eee eee reece (3 


for the drain pipe. The rate at which kinetic and potential energy is bei 
lost from the system, in the form of draining liquid is 


: w Ag x2 Ww 
- ——*— +— Ap L 
x( 2¢ g 2 1) 


and is equal to the rate of change with time of the kinetic and potential ener 
as given by Eqs. 33 through 36: 


dt (K.E.+P.E.) = a[% Ay (#- z,*) (a, ) x 


w AQ 2 Ag \2 Ware 
+ 26 Lx*+wAyj H- xx) a eee 
w Ag é 
= 2¢ +Ww Ag Lyx eee a ee ge ee er (3 


Eq. 37 yields 


Ag A2\* 1 Ag \2 A 
aches pad (et eae ROM IE esl PES ay ie “2 = 
[+52 (32) x Ju +9 [- (22) |s g [14 - x | =. 0.00 


Eq. 38 is identical to Eq. 5 with n H = F = 0, which was derived from mome! 
um considerations. 
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PROBLEMS CONCERNING USE OF LOW HEAD RADIAL GATES 


Discussion by Robert G. Cox and Arthur Toch 
Closure by Thomas J. Rhone 


OBERT G. cox,! M. ASCE.—The author is to be complimented for his 
lent paper covering the problems involved in the design and operation of 
ead radial gates.. The writer is particularly interested in the discussion 
scharge characteristics. 

1e fourth method presented by Mr. Rhone appears to be identical to the 
sdure, used by the U.S. Army Engineer Waterways Experiment Station, to 
lop a discharge coefficient chart for radial gates on curved spillway 
‘s. This chart was prepared for and included in “Hydraulic Design Cri- 
” sponsored by the Office, Chf. of Engineers, Dept. of the Army.“ The 
chart (Fig. 1 herewith) and Fig. 6 of the author’s paper are quite simi- 
2xcept for omission of the data points on Fig. 6. 

1e boundary geometry of flow under a radial gate on a curved crest is 
lex. To simplify discharge computations many engineers have used gate 
ngs and heads based on the gate seat and crest elevations. Furthermore, 
st instances, the effects of the direction of the approaching streamlines 
been ignored. Consequently, the results are often misleading when ap- 
to other structures of dissimilar geometry. 

1e WES method involves the net gate opening between the gate lip and the 
sd crest. It also includes the head to the center of the gate opening as- 
ng theoretical orifice discharge and atmospheric pressure under the is- 
‘jet. Discharge coefficients were computed and plotted as a function of 
ngle formed by the tangents to the gate-lip and to the curved crest at the 
ate opening. It is believed that this angle controls the direction of the 
mlines in the immediate vicinity of the gate opening. The results should 
fore be applicable to geometrically similar structures. 

1e WES chart was developed using model and prototype data. The data 
s shown are for spillway crests shaped to the form of the lower nappe of 
e jet over a sharp crested weir. The effects of gate opening-head ratios 
e discharge coefficients were studied. No well defined families of curves 
be drawn. Therefore, an average curve was accepted for design pur- 
. 

r, Rhone’s Table 1 indicates unusually good agreement between the meas- 
prototype and model discharges. The tabulated discharge computed by 
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DISCHARGE COEFFICIENT (C) 
WHERE: : 
Fi 
ORMULA “ Ge=NET GATE OPENING 
Q=CGoBV2gH B=GATE WIDTH 


‘a 
ER SeS SIU RSUERLT LAA CATAIN SRSOR SVAN PRBS CARAR PER SASRASe 


SYMBOL TEST GATE SEAT (X/Ha) 


MODEL DATA 
v CW 80I, AVERAGE 0.000 
4 CW 801, AVERAGE 0.167 
v WHITNEY DAM, AVERAGE 0.127 


PROTOTYPE DATA 


CENTER HILL, AVERAGE 
FORT GIBSON, AVERAGE 0.137 
WOODS RESERVOIR DAM 0.34 


@eao 
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0.70 0.75 


H=HEAD TO CENTER OF GATE OPENING 


FIG. 1.—DISCHARGE COEFFICIENTS FOR RADIAL GATES. < 
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surth method for the 3-ft gate opening is 4,880 cfs which is about 11% 
ter than the measured prototype discharge of 4,400 cfs. The writer, in an 
apt to discover the reason for this difference, computed a discharge of 
!cfs. The latter figure is within 2% of the measured discharge. A net 
opening of 2.90 ft, an effective head of 30.75 ft, and a 8 of 73.1° were de- 
ined graphically from available data.3 A discharge coefficient of 0.681 
Fig. 1 was used in the computation. Thewriter would be interested in 
ig comparative figures used by the author for his computation. Similar 
yutations for the 6-ft and 10-ft gate openings were in close agreement 
discharges computed by Mr. Rhone which are about 5% greater than the 
sured prototype discharges. This difference may result from the fact that 
on Ferry crest-shape, downstream from the gate seat, is considerably 
2r than the WES shape. 
‘the study of model data, it was noted that the discharge coefficients in- 
sed rapidly for gate openings less than 0.1 ft model. For similar proto- 
gate openings the coefficients were essentially constant or showed a 
t decrease. This difference between model and prototype data is attrib- 
to inherent difficulties in accurately measuring small gate openings and 
;in models. Therefore, only model data with gate openings greater than 
t were used in preparation of Fig. 1. 
he writer concurs with Mr. Rhone’s comments concerning the need for 
ting the effects of such factors as approach depth, pier shape, gate radi- 
runnion position, etc. However, it would also be necessary to determine 
ffects of numerous combinations of these variables. A less elaborate ap- 
ch for early usable results would be to standardize some of the geometri- 
actors in terms of the design head. 


RTHUR TOCH,‘ M. ASCE.—As noted previously,° it is believed that a 
ral mathematical analysis for the radial gate in a flat-bottom, rectangu- 
ection is not practical. Furthermore, it is felt that continued use of 
yn’s article (3) asa reference in the literature on Tainter gates is re- 
able. Horton not only chose to equate a contraction coefficient to a dis- 
ze coefficient but did so for gate shapes of quite dissimilar geometries. 
the four methods presented by Mr. Rhone, only the last appears to dif- 
rom the others. The first two formulas are evidently the same, as the 
ant value (2/3) V2g can be incorporated in the discharge coefficient 
ut any loss of generality for engineering applications. Both C’s must be 
sted, therefore, to vary in like manner, if the geometrical characteris- 
f the efflux boundaries are comparable. 

ie third formula differs from the first two only in the inclusion of the 
1/sin 9. It seems questionable to assume that inclusion of the lip angie 
s manner is of significance. Also, it is apparent that Fig. 5 merely re- 
the cosecant and cosine of any given angle. 

cause of the algebraic difference between the first three formulas and 
murth, the coefficient of discharge, for the latter, will have a quite dif- 
t numerical magnitude; it must still be dependent, however, on the same 
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variables or parameters. In fact, the importance of this aspect of the proble 
can hardly be overemphasized: The coefficient of discharge must be expect 
to be a function of the geometric configuration of gate and dam crest. In th 
respect, then, none of the methods presented by the author appears likely 
shed much light on the conditions of efflux or to increase the precision 
discharge prediction. 

In connection with Mr. Rhone’s request for further information regardi 
tests on radical gates, a study® made at the Iowa Institute of Hydraulic R 
search should be cited. Although this work did not include as many variabl 
as the author thinks necessary or desirable, Z. M. Glowiak did discuss a 
arrange the results he obtained ina rational manner. Another report® co 
tains fairly complete data on flow conditions and pressure distributions 0 
tained on a model of a radial gate on a relatively high dam. 


THOMAS J. RHONE, ! M. ASCE.—Messrs. Pariset and Michel have broug 
out some interesting facts in their discussion, and their extension of the c 
efficient-of-discharge curves is, indeed, a valuable addition. None of t 
crest profiles used in obtaining the curves presented in the writer’s pap 
was as thin as the Edgard de Sonza Dam crest profile; in fact, the highe 
free-flow coefficient was only 3.84. The larger values of the controlled fl 
coefficients, and the fact that for a given H the coefficient becomes essenti: 
ly constant when @ is greater than 70°, rather than increasing, as shown 
Fig. 6, indicates that other variables should possibly be considered in 
analysis. The use of the cavitation coefficient presented is a very good mea 
ure for potential cavitation damage. A similar equation8 has been used as 
means of determining cavitation potential in conduits and entrances, and 
the case of gate slots, has been proved with model prototype comparison: 
The writer assumes that the statement made by Messrs. Pariset and Mich 
“To determine the maximum allowable pressure on the profile of a spillway 
is in error, and should read either “minimum allowable pressure,” or “ma 
imum allowable subatmospheric pressure,” as stated later. 

Mr. Bowman’s awareness of the problems presented in the paper and |] 
subsequent comprehensive discussion are very gratifying. The increas 
number of papers that are appearing in technical publications seem to in« 
cate a favorable trend in “sharing the knowledge gained from basic researc! 
It is hoped that this policy will continue, and will be as well received as 
this instance. 

Mr. Bowman’s discussion of the various discharge quantity equations p1 
vides a clear, concise analysis of the advantages and disadvantages of ea 
He also provides a guide for the selection of the best combination of equatic 
to use under various conditions of gate design, gate opening, and head abc 
gate seal. His own method of obtaining an approximate rating curve is an : 
cellent “short-cut” system, and is probably as accurate as any of the me 
ods provided by the writer. The corrections for velocity of approach and ci 


6 “Deflusso Sopra Dighe Tracimate, Somontate de Paratoie a Settore,” (Flow 0 

Sees Surmounted by Sector Gates), by V. Calderini, L’Energia Elletrica, Milan, 1 

938. : 

7 Hydr. Research Engr., Div. of Engrg. Lab., Bur. of Reclamation, Denver, Colo. 

8 “Engineering Hydraulics,” edited by Hunter Rouse, John Wiley and Sons, p. 30. 

9 “Hydraulic Characteristics of Gate Slots,” by J. W. Ball, Journal of the Hyc 
Div., Proceedings, ASCE, Volume 85, No. HY 10, October, 1953, pp. 81-114. 
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tions, etc., are, for the most part, a matter of judgment, and will vary 
ding to the practice of each design office. 

‘he datum profiles used in the series of tests described in the paper were 
d on the hydraulic model investigations of the Hoover Dam spillway. 10 

Ir. Bowman’s discussion on gate seals is very informative. His observa- 
3 on prototype installations broaden the scope of the information provided 
1¢ paper, which, as he stated, was principally derived from Bureau of 
lamation practices. The writer cannot answer the question as to how the 
2al negotiates the bottom corners of a radial gate. In most of the struc- 
s investigated, the side seals were of the molded angle type, and the bot- 
seals were usually of the rectangular type, with a tightly fitted butt joint 
e corners. 

1 his discussion, Mr. Cox has brought forth an unfortunate omission; 
ely, that the writer did not identify the source of Fig. 6. The curves were 
ined from the Waterways Experiment Station chart, as he stated. The 
harge quantity that he questioned was given as 4,880 cfs in the table, but 
id have been 4,580 cfs, which is within 4% of the measured quantity. The 
er’s figures in arriving at this quantity were as follows: 


Net gate opening 2.96 ft 
Effective head 30.77 ft 
B CARES 
Discharge coefficient 0.681 


he small discrepancies between Mr. Cox’s values and the writer’s can 
ly be attributed to different interpretations of the graphical portions of 
solution, and either could be accepted. Mr. Cox’s explanation for the 
er values for the computed discharges at bigger gate openings is entirely 
ible. A steeper crest would permit higher flow quantities, other conditions 
g equal, thana flat crest. This is particularly true for the larger gate 
ings, as shown by Messrs. Pariset and Michel, and Bowman. The writer 
es with Mr. Cox’s statement concerning the standardization of some of 
reometrical design. The use of the design head as a basis for obtaining 
snsionless ratios is successfully used in other phases of spillway design, 
should be equally satisfactory in radial gate design. 

Ir. Toch implies that the discharge formulas presented in the paper are 
recommended by the writer; this was not the writer’s intention. The four 
10ds presented in the paper are inmore or less common use, and the pur- 
-in analyzing and discussing them was to show that each method is inade- 
e to some degree. The similarity between the first two formulas is in the 
of the equation only. It should be obvious that the coefficient-of-dis- 
ge values in the first equation are average values obtained from numer- 
model studies, irrespective of gate geometry. In the second equation, the 
ficients are related to the geometrical configuration of the gate. Appar- 
,, the originator of the third method believed that the use of any gate ona 
t would affect the coefficient of discharge to the extent shown in the lower 
hon Fig. 5, and that the gate geometry would further modify the coeffi- 
- in the manner shown in the upper graph. The writer cannot agree that 
pper graph “merely relates the cosecant and cosine of any given angle,” 
) “Studies of Crests for Overfall Dams,” Boulder Canyon Projects, Final Reports, 
Vi—Hydraulic Investigations, Bulletin 3. 


36 March, 1960 Hy 


as stated by Mr. Toch. It seems equally apparent that the graph presents 
simplified method for obtaining the value of 1/sin @ by using the read 
available values of the gate radius, trunnion height, and gate opening abe 
crest. Mr. Toch is correct in believing that the coefficient of discharge 
also a function of the crest shape. Previous discussions have vividly demo 
strated this to be true, and it should be an important consideration in ge 
eralizing a radial gate discharge equation. The two references listed by N 
Toch are appreciated. The thesis by Mr. Glowiak has recently been review 
by the writer, and it is a competent piece of research work that adds cc 
siderably to the knowledge of radial gate discharge characteristics. 
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DISCUSSION OF RELATIONSHIPS BETWEEN 
PIPE RESISTANCE FORMULAS 


Closure by Walter L. Moore 


VALTER L. MOORE, |! M. ASCE.—Several of the discussors, recognizing 
inconvenience of the Colebrook-White relations for solving resistance 
lems, presented different approaches to aid in the achievement of a more 
renient method. Emphasis was given to the limitations of the Colebrook- 
fe transition function for unusual roughness forms and thus contributed to 
ore comprehensive presentation of the resistance problem. Some of the 
ements in the discussions indicated a misunderstanding of the paper and 
e will be cleared up to the extent possible. 

fr. Ackers’ discussion was particularly interesting in that he described a 
rent attempt to present the Colebrook-White relations ina more usable 
a than previously available. The first of the two papers to which he refers 
s an excellent comprehensive review of the pipe-resistance problem. In 
mce, the approach described by Mr. Ackers consists of plotting the Cole- 
k-White relations in terms of dimensionless velocity, size, and slope 
meters. It is of interest tosee how the parameters he used can be formed 
the combination of the more conventional dimensionless parameters of 
tion factor f, Reynolds Number Np, and relative roughness e/D. Ackers’ 
snsionless parameters are listed below, together with their definition from 
oaper and their equivalent in terms of conventional parameters: 


V=Ve/v = Nr 
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ckers’ Eq. 1 the k is apparently equivalent to e as used in the rest of the 
iT. 

‘he general resistance diagram, plotted in terms of the foregoing para- 
ars, represents a logical arrangement and is based on sound dimensional 
oning. It has certain advantages which are readily apparent for problems 
lving a simple single pipe. A direct solution may be obtained for either 
yelocity, the diameter, or the head loss, if the other necessary quantities 
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are known. It is, thus, definitely more convenient than the Moody diagram f 
problems in which the diameter, the velocity, or the discharge are unknown 

The new diagram proposed by Mr. Ackers is not an improvement in all r 
spects however. Since each of the parameters V, R, and S varies over a wi 
range of practical values all of the scales must be rather compressed, th 
impairing the precision with which values can be read from the charts. 
course, large charts can be made to improve the precision (see charts int 
pouch of Ackers’ reference 7) but these become rather awkward to use a 
are not likely to be widely available. On the other hand, the pipe-resistan 
factor f varies over a rather narrow range of numerical values, and in t 
Moody diagram the scale for f can be chosen large enough to read with t 
necessary precision. 

The physical significance of the variables is not as readily apparent int 
Ackers’ diagram as it is in the Moody diagram. The interrelationship of t 
velocity, size, and head loss (or slope) is apparent for the proposed diagra; 
but this is already apparent from the general form of the various resistan 
relations. The significance of changes in the roughness or the fluid viscos 
are not readily apparent from the new diagram, because these variables o 
cur in more than one of the dimensionless parameters. The Moody diagra: 
on the other hand, clarifies the significance of these variables and brings ¢ 
more Clearly the meaning of the “smooth pipe” and “ rough pipe” limits. 

For complex pipe systems, including transition losses or pipes in serie 
the writer can see little advantage in the Ackers diagram as compared w 
the Moody diagram. If the rate of discharge or the velocity is known, eith 
diagram may be easily used. If the discharge is not known a trial solution w 
be required with either of the diagrams. 

It is for these complex pipe systems, as well as for pipe networks, that t 
writer has found the exponential relations so useful. The writer believes it 
more rational to examine the range of Reynolds Number and relative rou 
ness likely to be encountered and to develop a suitable exponential resistan 
equation from the modified Moody diagram, according to the methods he f 
proposed, than to arbitrarily “select” an exponential formula without due co 
sideration for the range of parameters involved in the particular problem. | 

Mr. Szesztay briefly reviews his work applying the Colebrook-Wh 
formula to a graphic solution of channel flow. Essentially his method consi: 
of specifying the shape parameters of the cross section in terms of the wat 
depth h, thus making the area and hydraulic radius a function of water deptl 
only. For a given sand-grain roughness k the discharge may then be « 
pressed as a function of water depth h and slope by means of the Cru 
modification of the Colebrook-White formula, as given in Szesztay’s Eqs 
and 4, The chart presented in Fig. 3 permits a ready solution of these rel 
tions. One could easilyalter the relations given by Mr. Szesztay for the cha 
nel shape as a function of water depth, if he wished to use different crite: 
for determining the channel shape. 

Although Fig. 3 is clear for the assumed k-value, it is not clear to 
writer how the corresponding value of n = 0.025 can be stated. It would se 
that the value of n corresponding to agiven roughness value k would bea fur 
tion of the hydraulic radius and, hence, of the water depth h. Mr. Szeszi 
states that statistical studies of n and k apparently fail to reveal any c¢ 
sistent effect of the value of the hydraulic radius R. It would seem that t 
result would be worthy of additional study. 
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in Fig. 3, the key for the slope S is listed as centimeters over centimeters 
‘reas it appears that it should be centimeters over kilometers. 

fhe writer was pleased to note the similarity of viewpoint presented by 
Miles. His extension of the concepts to a tabular comparison of e/D and 
en and Williams’ C-values was particularly valuable. The results that the 
en-Williams C approaches 160 as the relative roughness approaches 0 was 
resting. The reasonableness of this value is indicated by the C of 155 for 
bituminous enameled-lined pipe of 16 in. or larger given in the 5th edition 
3abbitt and Doland’s “Water Supply Engineering.” 

The writer is grateful to Mr. Bilonok and Mr. McPherson for correcting a 
1erical error in the first example of the appendix. The figures shown give 
ilue for hg of 4.88 ft rather than 5.84 ft, as printed inthe paper. The cor- 
t figure results in a value of Kp, = 0.455 thus giving a resistance equation 
this example of: 


he = 0.455 Q1-88 


3 numerical erroralso led to an incorrect value for the lengthof the 12-in. 
1 pipe appearing in Table 2. In order that the equation for the 12-in. steel 
, he,» = 1.99 Q1.88 be correct, the iength of the pipe should be 4,370 ft in 


-e of the 3,660 ft shown in Table 2. 

Mir. Bilonok also recognized that exponential formulas are in common use 
yracticing engineers because of their convenience, as compared with the 
eral resistance diagram and the Darcy-Weisbach formula. 

n discussing the effect of aging on pipe resistance, Mr. Bilonok states that 
shness factor e, is a “complete measure of the size, shape, and distribu- 
of individual roughness elements.” Although the introduction of roughness 
ameter e indicates a great improvement over other methods of charac- 
zing surface roughness, the writer agrees with Mr. McPherson and Mr. 
ers, that shape and distribution of individual roughness elements are not 
assarily characterized by the value of e. For the case of statistically dis- 
uted roughness, the single parameter e appears to be a satisfactory char- 
ization of the entire roughness. It has been shown that in the case of ex- 
ne variations of roughness form, such as grooved and corrugated pipes, 
e other geometric parameters are important besides a measure of the 
shness heights. 

"he writer is puzzled by Mr. Bilonok’s question regarding the finding of 
exponent in example 2 where m = 1.95, which he says, “is not in agree- 
t with fand e/D.” In the second example in the appendix if the values of 
and f, or e/D listed are used to enter Fig. 1 of the paper, the value of m 
be found to lie between 1.95 and 1.96. 

‘he writer is indebted to Mr. McPherson for a very painstaking and com- | 
1ensive discussion which stressed some points which probably deserved 
e emphasis. Also, Mr. McPherson apparently misunderstood a number of 
ements in the paper. ; 
tis obviously true that no system for approximating a complex relation- 
can be more precise than the relationship being approximated. Thus the 
\osed exponential equations are certainly no more precise (but fortunately 
slightly less precise within the limitations imposed) than the Colebrook- 
-e relationship, from which they were derived. Mr. McPherson states that 
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“the preciseness of the general resistance diagram implied by the ‘exponen- 
tial equations’ of the author does not appear to be justified.” It is a misunder- 
standing of the original paper to consider that the exponential equations imply 
about the preciseness of the general resistance diagram. They are merely an 
attempt to put the information in the general resistance diagram in a more 
convenient form for a number of types of resistance problems. 

Mr. McPherson is concerned because the numerical values in Table 2 pass 
beyond the realm of water-works design practicality. 

The table was designed merely as a means of illustrating the effect of 2 
difference in the exponent m on the concept of pipe equivalents. Of course, the 
values extend beyond the range of water-works practice. It was not the 
writer’s intention that the concept be limited only to water-works practice. 

Mr. McPherson questions the statement that the exponential equation will 
agree with the general resistance diagram within 5% for a 20-fold variation ir 
Q. This statement is significant, not because one expects to analyze a wate 
distribution system over a 20-fold variation in Q, but in that it indicates that 
a satisfactory exponential equation may be developed even when the initially 
expected value of Q is poorly chosen. It must, of course, be realized that fo1 
the exponential approximation to apply for a full 20-fold variation, it must be 
developed about a central point in the range rather than about one extreme 
end of the range. Thus, in the first example in the appendix, which gives ¢ 
head loss of 4.88 ft at a Q of 3.53 cfs, the comparison should be made at 


and at Q = 3.53 V20 = 15.8 cfs. 
From the equation hg = 0.455 Q1.88 
he = 0.290 ft at Q = 0.79 cfs 
and hg = 82 ftatQ = 15.8 cfs. 
From the general resistance diagram 
hg = 0.299 ft at Q = 0.79 cfs (f = 0.019) 
and hg = 85 ftat Q = 15.8 cfs. (f = 0.0135) 
The percentage differences are as follows: 
at Q = 0.79 cfs; (0.292/0.299 - 1) = -2.2% 
at Q = 15.8 cfs; (82/85 - 1) = -3.5% 


The writer agrees with Mr. McPherson that in designs of water-works net 
works a constant value of the exponent m may normally be used. The write 
believes, however, that an exponent chosen from the modified Moody diagrat 
is more rational and will also give better results than arbitrarily using th 
1.85 exponent from the Hazen-Williams formula. Mr. McPherson’s discussio 
does not seem to recognize that network problems also exist in fields othe 
than water-works practice. In some instances it may even be desirable to us 
different exponents m for different branches of the system. Contrary to M 
McPherson’s statement, digital computers are not necessarily restricted 
using exponent m.2; 
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2 “More on Pipe Line Computers,” by Robert L. McIntire. Oil and Gas Journal, V% 
57, No. 26, June 22, 1959. 
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Mr. McPherson has apparently not understood the writer’s use of the term, 
‘exponential resistance equations,” which was defined in the paragraph im- 
nediately following Eq. 2. Thus, any equation of the form hg = K Q™ (which 
ncludes the Hazen-Williams formula as a special case) is considered an ex- 
onential formula. Mr. McPherson must have misunderstood the term in 
irder to make the statement that he was “forced to reject as obviously un- 
‘ealistic (and clearly not proven) the statement by the author: ‘although the 
nethods for analysis of networks will not be considered here, it is pertinent 
o note that exponential resistance relations will be most convenient for use in 
his type of analysis.’” Mr. McPherson goes on to propose the use of the 
Jazen-Williams formula for network analysis which is really a confirmation 
f the writer’s statement although restricting it to the special case of an ex- 
ponent m = 1.85. 

Again Mr. McPherson apparently misunderstood the statement that, “thus 
he Hazen-Williams formula with its exponent of 1.85 corresponds to a com- 
ination of relative roughness and Reynolds number lying along the line, 
n = 1.85 on the modified Moody diagram.” The discussion preceding this 
tatement clearly deals with the effect of m on the change of head loss with 
hange in Q. The statement should thus be taken to indicate that only along the 
ine m = 1.85 will the Hazen-Williams equation have the correct exponent. 
obviously the whole field of friction factors and Reynolds numbers may be 
overed by lines corresponding to an exponent m = 1.85. This is what is done 
n the chart presented by Mr. McPherson as Fig. D-2. The important point is 
hat the change in f with Reynolds number occurs at the same rate for both 
ne Hazen-Williams equation and the general resistance diagram, only along 
ne line corresponding to m = 1.85. This may be clearly shown by plotting 
oordinates of the line m = 1.85 from the modified Moody diagram onto Fig. 
\-2. It will be seen that only along this line are the solid lines representing 
le Hazen-Williams equation approximately parallel to the dashed lines 
epresenting constant relative roughness on the general resistance diagram. 
'o the right of the chart the Hazen-Williams equation would indicate a con- 
nuing decrease of friction factor with increasing Reynolds number, whereas 
1e rough pipe relations represented on the resistance diagram show that the 
ciction factor should become constant. It might be argued that Cy yw, could 
e decreased as Reynolds number is increased, thus forcing Hazen-Williams 
uation to yield a constant f at high Reynolds number. This would complete- 
r destroy the added convenience of the exponential type equation in this 
ange. Clearly what is needed in the rough-pipe range is an equation with an 
sponent m = 2 in order to fit the conditions for a constant f-value. 

Further evidence of the correctness of the original statement is given by 
e first equation presented by Mr. McPherson under the heading, “The 
azen-Williams Formula.” This equation 


: 1060 

HW. = — = 

0.15 70.015 
Cyw. Np * D 


the same form as Eq. 5 of the paper with m = 1.85 because D0-015 will be 
actically 1 for any reasonable value of D. Thus if Co, is set equal to 
160/ Cy w, and m = 1.85, Eq. 5 will be seen to correspond to the foregoing 


uation. | 
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The data presented by Mr. McPherson in Tables D-3 and D-4 for nev 
linings and old linings is of interest. If the C-values given for Reynolds num: 
ber and friction factor are plotted on the modified diagram, the points for the 
new linings, for the most part, fall between the lines of m = 1.90 and m = 1.96 
For the old linings (Table D-4) the points, for the most part, lie between th 
lines-m = 1.94 and m = 1.975. This data, also, indicates that the resistance 
relations for the older pipes would be better approximated by a larger valu 
of the exponent m than for the new pipes. It is also of interest that the indi. 
cated values of m are all considerably about the 1.85 line corresponding t 
the Hazen-Williams equation. Thus this data also supports the writer’s con 
tention that an exponential equation with an exponent larger than 1.85, in fac 
one approaching 2, would be preferable for old water mains. 

It is the writer’s belief that the modified Moody diagram in conjunctio 
with available data for e, or better still test results such as presented by Mr 
McPherson, will enable the engineer to select the exponential equation whic! 
will be most appropriate for the type of problem with which he is dealing 
The writer recognizes that there is a possibility (which at this time seem: 
rather remote) that large changes in roughness form due to encrustation an 
the forming of ridges or waves might alter the form of the transition functio: 
from “smooth” to “ rough” pipe conditions and thus affect the proper value o 
the exponent m. That the shape of the roughness should change enough to sig 
nificantly affect the treatment function seems unlikely, and in the absence o 
reliable experimental data on this point, the modified Moody diagram appear 
to be a reasonable improvement over existing methods. 

Additional research to explore the shape of the transition function and thi 
geometric form of roughness for encrusted pipes would certainly be of inter 
est. In such an investigation it may well be that the reduction in diameter wil 
be more significant than the change in roughness since the diameter enter 
into the equations to the fifth power. 

The work described by Dawsey attempting to relate measurements of rough 
ness profiles to the equivalent sand grain roughness e has much appeal. T 
be able to run a test on a small portion of a conduit lining, or better still on 
cast taken of a portion of the lining, and convert this to a hydraulic roughnes 
parameter would be a tremendous forward step. There are many difficul 
problems of measurement, and meaningful analysis of the measurements, tha 
will have to be solved before this approach is likely to yield reliable results 
The magnitude of the potential gains are so great that efforts toward this goa 
seem well directed. The writer thanks Mr. Dawsey for noting the omission ¢ 
parentheses in the Colebrook-White equation as presented in Table 1 of th 
paper. 

The writer is grateful to all of the discussors for contributing addition 
information and new thoughts regarding pipe-resistance relations. The dis 
cussions have supplemented the paper by presenting other attempts to de 
velop convenient as well as reliable methods of evaluating pipe-resistance. | 
is hoped that the concepts presented will lead to an improved understanding 
pipe resistance relations. Time and the experience of the practicing enginee 
will determine which of the various approaches will prove to be most useful. 
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CHANNEL-SLOPE FACTOR IN FLOOD- FREQUENCY ANALYSIS2 


Closure by Manuel A. Benson 


MANUEL A. BENSON,! M. ASCE.—Mr. Sammons is to be commended for 
is thorough analysis and reanalysis of the writer’s data. The original article 
as intended to cover only thechoice of a simple, yet efficient, index of main- 
hannel slope. Mr. Sammons covered a wide scope, including basic data col- 
ction, mapping standards, flood-frequency analysis, multiple correlation, 
atistics, and the use of automatic computers. Much of his discussion con- 
erns points which will be explained in detail in the pending final report of the 
ympleted project. Mr. Sammons makes some sweeping indictments of vari- 
Is practices, to most of which the writer wishes to plead “not guilty.” For 
‘ample, the writer sincerely believes he does not share that misconception 
the use of automatic computers, and he is certain that he does not apply 
ndiscriminate weighting coefficients” to variables. 

Mr. Sammons puts much stress on the distribution of the variables which 
re used in the regression analysis. Multiple-correlation techniques, to be 
lid, require linear relationships between the variables. The transformation 
logarithms of the original variables was made on the basis of comparative 
garithmic and rectangular plotting between the dependent and independent 
\riables, and plottings of the residual errors at various stages against the 
dependent variables. One variable (surface-storage area in ponds) used in 
e final analysis was adjusted by a “ rectification constant” determined graph- 
ally to be necessary to linearize the logarithms of the variable. Finally (as 
iggested by Mr. Sammons) the cumulative distributions of the variables 
entually used were plotted on log-probability paper and were proved to be 
sentially linear, within 2 standard deviations of their means. This con- 
rmed the validity of the logarithmic transformations. 

No logical reason is apparent for Mr. Sammons’ statement that the curvi- 
lear regression constant of “a” in Fig. 3, indicate that “a curvilinear 
zarithmic model would best fit the data.” Elsewhere, Mr. Sammons refers 
the “concave upward bias of the constant term a on Fig. 3.” A bias is much 
ore apt to be introduced by an a priori theory than by adherence to the data. 

The distribution of most of the variables, which may influence flood peaks, 
sufficiently close to log-normal, thus that distribution may be assumed. In 
y event, their exact distribution is not critical for the purpose of relating 
2m to flood peaks. However, the writer is not willing to assume log-normal 
stribution for the annual peak discharges at a station, because this is a 
itical assumption. The form of the distribution becomes the principal basis 
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for determining extreme flood peaks. Statistical methods may be nonpara- 
metric as well as parametric. The writer prefers to correlate topographic 
and climatic variables with the “percentiles” of the flood-peak distributior 
rather than with its parameters (mean, standard deviation, and skew), simpl} 
because the use of parameters means making assumptions that are no 
necessarily or generally true and that may actually be in error. 

Mr. Sammons cannot suppose that the investigation which he cites, it 
which precipitation proved to be more important than the drainage area, i: 
generally applicable or reflects on the results for New England. If the analy: 
sis were made for a group of small drainage areas (necessarily limited i 
range of sizes) within a region of widely varying precipitation, it is to be ex: 
pected that precipitation will prove more important. In fact, if all the drain. 
age areas are the same size, drainage area would be of no importance what: 
soever, at least as shown by the multiple-regression analysis. 

The writer agrees with Mr. Sammons that the median values of variable: 
might be better indices than the means. However, the U.S. Weather Bureat 
publishes only the means for climatic factors such as precipitation and tem. 
perature. The mean of topographic variables, such as elevation and lan 
slope, is more simple to compute than the median. 

The scales of topographic maps of New England were adequate to defin 
profile elevations for the basins used in the analysis. Regions of flat slopes 
less adequately mapped regions, or extremely small watersheds might re 
quire surveys or photogrammetry to determine the channel slope. The write 
did not investigate the effect of variation in map scales because the availabl 
scales introduced no problems for the data studied. 

The ultimate findings have related flood peaks of recurrence intervals 
ranging from 1.2 to 300 yr, to the following hydrologic variables: 


. Drainage area; 

. main-channel slope; 

. surface storage in lakes and ponds; 

. precipitation intensity; 

. mean January temperature (an index of the effect of snowmelt on peaks) 
and 

6. orographic factor. 


oP WON 


A report on the investigation and the findings is nearing completion. 
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ROLL WAVES AND SLUG FLOWS IN INCLINED OPEN CHANNELS2 


Discussion by Tojiro Ishihara, Yuichi Iwagaki and Yoshiaki Iwasa 


TOJIRO ISHIHARA, ! YUICHI IWAGAKI,2 M. ASCE and YOSHIAKI 
WASA.3—The writers have reviewed this paper with interest, because we 
Iso have made continuous investigations of the hydraulic characteristics of 
oll waves and other associated problems in laminar and turbulent flows from 
heoretical and experimental approaches during the past several years. Un- 
ortunately, however, our results are quite different from the author’s study, 
nd we also do not follow his treatment of roll waves and slug flows from the 
nowledge we obtained through research in the hydraulics laboratory of Kyoto 
niversity. 

Our research project to disclose the hydraulic characteristics of roll 
aves, needed by the engineering request for soil conservation, have been 
tudied by the writers and others since 1950. Several reports were already 
ublished in Japanese4 -7 and in English8,9 after completing specific problems 
elated to roll waves. Through our research works, the basic concept of the 
‘eatment for the hydraulic characteristics of roll waves is essentially due to 
1e studies of H. Thomas10 and R.F. Dressler.11 More precisely, the deriva- 
on of the hydraulic characteristics of roll waves is based on the fact that 
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(1) the wave velocity of roll waves is definitely constant (terminal velocity by 
the author), regardless of their flow regime, when the final shape in wave 
pattern is obtained, and (2) the flow velocity is, in all cases, less than the 
wave velocity of roll waves in magnitude, thus the flow regime of rear side of 
a roll wave is shooting, whereas that of the front side, including the wave 
crest, is tranquil when the train of roll waves is observed from the moving 
coordinate system travelling at aconstant speed of wave velocity. The former 
condition was also verified by the author, while the latter condition of Thomas 
was not equivalent to the author’s in all points because the author divided the 
flow behaviors into two different wave patterns, depending on the value of the 
Reynolds number. 

In order to ascertain our standpoint concerning the author’s derivation for 
the hydraulic characteristics of roll waves and slug flows, our analytical 
treatment on the hydraulic characteristics of roll waves and the experimental 
verification by the direct measurement of roll waves by means of electronic 
devices will be first described: 

Based on the principle of a constant wave velocity indicated in the fore- 
going, the unsteady behaviors of roll waves can be reduced to those in steady 
regime by eliminating the time derivatives in the basic quasi-linear partica. 
differential equations of hyperbolic type for unsteady flows, when viewing the 
wave pattern from the moving origin. The hydraulic characteristics of rol. 
waves are then easily derived in the following. 

The equation of momentum for unsteady flows in open channels with a con- 
stant inclination is 


ov ov oD VidDirwed g v2 | 

Br to Ving + £008 8 ay + (la) ar 8 ein ee 
and the equation of continuity is 
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in which the pressure distribution is assumed hydrostatic, x is the distanc 
from the origin along the channel bed, t represents time, V is the meai 
velocity of flow, D denotes the depth of water, A is the flow area, R represent: 
the hydraulic radius, C is Chézy’s coefficient (assumed variable), g denote; 
the acceleration of gravity, 6 is the inclination angle of channel bed, and « 
equals the momentum correction factor usually defined. Apparently, th 
velocity profile of laminar flow is of the parabolic type which was already de 
rived by many hydraulic engineers as well as the author and the momentul 
correction factor becomes then aconstant of 1.2 with a simple calculation. O 
the other hand, the evaluation of a in turbulent flow cannot be made in a defi 
nite number, because it is still a function of flow and channel characteristics 
but it is satisfactorily assumed 1.05 in engineering practice. In the derivatio 
of Eq. 1, the curvature of a stream line is neglected, and thus the influence o 
surface tensions is also ignored, which is different from the author’s treat 
ment. Of course, the surface tension will have to be considered for lamina 
flow in very steep channels. : 

Taking account of the preceding principle for a constant wave velocity, th 
time derivatives in Eqs. 1 and 2 can be eliminated by setting =x - Vy 1 
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i the resulting equations obtained by transforming the independent variables 
mm x and t to ¢ are 


4 g(V - Vw)(dA/dD}{sin @ - (V2/Cc2 R)} 


& ~ (@A7dDX(a V - VwXV - Vy) +(1-a)V Vy} -gAcoso''*** (3) 


1 


g A{sin 9 - (v2/c2 R)} 


- (dA/dD}{ (a V - Vw)(V - Vw) + (1-a) VVy}-¢gAcos@°°*** (4) 


reality, the wave velocity may not be considered constant in a channel of 
‘egular section, but the treatment of analysis will be proceeded from a 
croscopic aspect for mean flows. Eliminating V from Eqs. 3 and 4, the 
adied wave pattern observed from the moving origin can be expressed by 


: g A(sin 6 -{(Vy A - K?/C2 R At}) f,(D) 


which K is a progressive discharge rate defined by (Vy - V)A =a constant. 
seen in Eq. 5, the wave pattern of a roll wave cannot be derived in a peri- 
c form observed in laboratories, and thus the solution for roll waves must 
obtained as a discontinuous wave pattern by combining each of water sur- 
€ profiles through the shock condition as Dressler did. Dressler also ob- 
ned a continuous periodic solution of roll waves being in forms of the 
iptic function as the second order solution of the shallow water wave theory. 
e actual wave pattern, however, observed in channels is rather in a similar 
m to discontinuous periodic solutions. 

The basic concept of Thomas obtained by his specified experimental works 
icates that the flow involves a control section in each water surface pro- 
1, as the flow regime changes from tranquil to shooting, so that in Eq. 5, 
h numerator and denominator become simultaneously zero and all the hy- 
wlic characteristics of roll waves can be uniquely determined at this sec- 
1. The resulting expressions for the mean velocity and the hydraulic radius 
the water depth in specified cases at the control section are 


Vo a(dA/dD), - ¥ a(a - 1)(dA/dD)2 +{ (dA/aD).8,/Net 


we en (6) 

Vw o(dA/AD)g - (So/N2,) 
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gcos@ 1 | a(dA/dD)g - Va(a- 1)(dA/dD)6 +4 (dA/dD) 9So/NFt |“. - (7) 
_ EE ed ee 1 SS as 
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vhich S is the (ey perimeter, Np denotes the Froude number defined by 
(g Ro cos 9) 2, and the subscript o indicates the values at the control 
tion which are eventually equal to the quantities in uniform flow. The pro- 
ssive discharge rate can also be calculated by the use of Eqs. 6 and 7. 
ff the flow is laminar in a rectangular channel, the normal mean velocity 
low, the normal depth of water, and the progressive discharge rate are 
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explicitely indicated as functions of the wave velocity and Froude number |: 
follows: 


Vo (6/5) - ¥(6/25) + (1/N2,) 
Vw (6/5) - (1/N2, ) 


(6/5) - ¥ (6/25) + (1/N2) : 


Do g cos 6 1 


ya We, (6/5) - (1/N2.) 
and 
K gcos@ 1 [ (6/5) - ¥ (6/25) + (1/2) }2 
= ah 
y= Ne (6/5) - (1/N2) 


ae - V (6/25) + (1/N2,) 


(6/5) - (1/N2.) 


Water-surface profiles as each wave pattern of roll waves are obtained | 
integrating Eq. 5 upstream and downstream from the control section. A 
though water-surface-profile equations cannot be integrated analytically 
general cases, those for laminar flows and for turbulent flows, characteriz: 
by a constant Chézy coefficient in rectangular channels, are integrable. 
cases of two-dimensional laminar flows, the resulting water-surface profi 
is indicated in an inverse functional form as follows: 


1 D2 cos a+ (6/5)(K?D,/g D2) + (6/5)(K?/g Do) eens iD 
sin 0 DA - Dp Sed - Da 


D2 cos 6 + (6/5)(K2Dp/g D2) + (6/5)(K2/g Do) merece: 
PatiPR if: Do - Dp 


+ cos 6(D - Do) 


in which the integral constant is selected so that D = Do at £ =0, and Da a 
Dp are 


and Do > Da > 0 > Dg. 
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The complete discontinuous solution of roll waves can be obtained by com- 


ing each of water-surface profiles through the shock condition, as will be 
scribed in the following. 


The shock condition for open-channel flows involving a moving discontinu- 


is expressed by the following three conditions: 
Mass conservation: 


-p Ap(Vy - Vp) = - p Ag(Vy - Vp) = Mo... eo... (14) 


Momentum conservation: 


Vw - V») ~ (Vy - Vp) =pgcos 6{ Ag(De - Dep) - Ap(Dp - Dep! ee er (15) 
Energy dissipation rate: 

_qdE MOM gcos 6;,2 

ON 2A; Ap td (Dp = Dgp) - Ap Az (Dp + Dgp) 


+ Ap Ag (Dp + Doe) - AZ (De - Day}... 2 (16) 


Thich Dg is the distance from the channel bed to the centroid of flow area, 

the subscripts of f and b indicate the values at the lowest and highest 
its of water-surface profiles. Using the expression for progressive dis- 
rge rate and Eqs. 14 and 15, the relationship between the lowest and high- 
depths of a wave becomes 


K{(1/A,) - (1/Ag} = g cos 6{Ag (Dz - Dgy) - Ap(Dp - Dgp)} ...- (17) 


oting the wave length of a roll wave by L, the following relationship must 
it between two successive waves: 


21D Sy SD) ee ee (18) 


at the shock front 
5 Dy(Dp) = 5 Day 1 (Pf) Tent oe et de eek a (19) 


[ the water-surface profile for a particular roll wave is obtained by inte- 
ing numerically Eq. 5, the complete solution of discontinuous periodic 
2tn of roll waves will also be calculated graphically or numerically. In 
‘Ss of two-dimensional laminar flows, the wave length of a roll wave is ex- 
tly determined by 


Dj cos 6 + (6/5)(K2Da/g D2) + (6/5)(K2/g Do) og, Db_DA 
Da - Dp Be De - Da 


1 
sin 0 


D2, cos 4 + (6/5)(K*Dp/g D3) + (6/5)(K?/g Do) Pages Dp 
Da - Dp Se D; - Dp 


+ cos 6(Dp - Dr) 
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The criterion for roll waves to maintain their final wave pattern, as d 
scribed in the foregoing, will be studied next. It is of common observati 
that the wave profile of a roll wave is monoclinally increasing from the re 
side to the front side, so that the following relationship for the water dey 
must be obtained at the control section, when the waves are observed fr¢ 
the moving origin. 


With the aid of Eq. 5, the above condition for laminar flows characterized 
the parabolic velocity profile and a = 1.2 becomes 


or, expressing in terms of Froude number, 
oc (=) Sh : | 
NE 2 So/\ ap ee m* = 0.8 m = 0.2) «2.3. .55 sie ( 


in which m is the shape parameter defined by 1 - Ro(ds/dA)o and it gradua 
increases to unity with increase in width of a channel. 

If the flow is turbulent, the criterion is also calculated in the same mz 
ner, and it is 


m Vo (2 Ry, /dC 2 
ee 2 (— } +— +1}2 1 Te ee ( 
2Vy (Co \dR/g m 


for the channels in which the flow area is proportional to a power of the wa 
depth. To obtain the criterion in a more familiar fashion for hydraulic en 
neers, the velocity formula of Vedernikov’s type will be introduced as 


V8 op Rd*? ein O75.) oe ( 


in which a, b and f are constants determined by given flow and channel chi 
acteristics, so that the Chezy’s coefficient C becomes 


oe pl/a Ry2(1+b)-a/2a (sin 9)2-a/2a 


Inserting Eq. 26 into Eq. 24, the resulting condition becomes 


or, expressing in terms of Froude number, 
dA 1+b)\2 | 
Nz, 2 So/ (SB) jm (++*) - 2(a- 1)m (+4) - (a- »| oe 


Eq. 27 is evidently the Vedernikov number for the initial instability of she 
ing flows. The last of the writers? also derived the same criterion by dec 
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the initial condition for continuous time growth of an infinitesimally small 
turbance with the use of Eqs. 1 and 2. A. Crayal2 and Dressler13 studied 
) similar criteria, and in specified cases of particular velocity formulas 
freysl4 and Keulegan and Patterson!5 obtained the instability conditions. 
most significance, in deduction of the foregoing maintenance conditions for 
l waves being in their final wave pattern, is that two conditions for main- 
ance of roll waves and for initial instability of shooting flows are exactly 
/ Same expression of Vedernikov, regardless of essentially different ap- 
aches in deduction. The formation of roll waves, therefore, will result 
m the instability of open-channel flows as many hydraulic engineers de- 
tibed. 

The previous explanation for the hydraulic characteristics of roll waves is 
ummary of our past analytical treatment developed by means of the basic 
cepts of Thomas and Dressler. The experimental works for the charac- 
istics of roll waves, to clarify the hydraulic behaviors of waves in definite 
ms and to verify the theoretical deduction of roll wave characteristics, 
re also paralleled by the writers in the hydraulics laboratory of Kyoto 
iversity. Two test flumes of a rectangular section were used for this pur- 
3e, and in the earlier stage of research progress, which was mainly sub- 
t to laminar roll waves, the aluminum channel of 20 ft in length was used, 
ereas the experimental runs for turbulent roll waves, which are calledas 
g flows by the author, were made in the lucite flume of 36 ft in length and 
variable slopes from horizontal to 30°. Continuous recordings of wave pro- 
1s of roll waves and of wave velocities were made by the electronic wave 
orders of resistance type, and the discharge measurement was made by 
ume metering of water at the downstream end. The details of experimental 
yeedures were reported in the writers’ paper.4 Comparing the theoretical 
sults described in the foregoing, the experimental behaviors of roll waves 
ained by the writers will be presented in the following. Data of the author 
| Binnie16 will be also supplemented. 

The criterion for initial instability and formation of roll waves will be 
ated first. For the sake of simplicity in illustration, the flow is assumed 
-dimensional. In this case, Eq. 23, which indicates the criterion for lami- 
‘flows, becomes 


NRA Leigh ky ITAL A (29) 


which Ry, is the Reynolds number given by Vo ‘Daf v, and the same ¢ri- 
ion for turbulent flows is also in terms of Froude number, 


nz, = 1/ { (222)° - 20-9 =) -(a- 0} beset (30) 
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and consequently, it depends on the velocity formulas like Chézy’s or Ma 
ning’s and the momentum correction factor. Since the Froude number is « 
pressed as a function of the Reynolds number, Manning’s roughness coef 
cient, slope and kinematic viscosity, when the Manning formula is used, 1 
instability criterion can be graphically illustrated in terms of the slope a 
Reynolds number as parametric expressions of the roughness and the m 
mentum correction factor. Fig. 1 illustrates the criterion for instability 
flows and formation of roll waves obtained in the preceding manner under t 
assumption that the kinematic viscosity is 0.01 sq cm per sec. Two zones 
stable and unstable flows will be apparently observed in the figure. In ope 
channel flow, the transition from laminar to turbulent flows takes place 1 
Reynolds number from 500 to 1,500, as often indicated by many enginee1 


= Laminar Transition ed Turbulent G:. 
peau O O ol e , 
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FIG. 1.-REPRESENTATION OF CRITERION FOR FORMATION OF ROLL WAVES. 


The second of the writers’ found that the Manning formula was not suitable 
the transition region, but (1 + b)/a was commonly very small and it becar 
negative in particular cases depending on the channel slope, through the e 
perimental research. The formation of roll waves, therefore, will not resul 

Now, let us consider the behavior of roll waves in formation when the fl 
is carried inan openchannel. For clarification in explanation, the chan1 
slope is assumed 0.02. When the regulating valve for discharge is open 
from closure, the flow begins to be carried. The roll waves are not observ 
until the Reynolds number becomes 50, as previously indicated. With i 
crease in discharge, the laminar roll waves become appreciable. Once t 
Reynolds number increases to 500, so the transition from laminar to turk 
lent flows begins to take placeand the roll waves disappear. This phenomen 
may be corresponding to that called the confused pattern by the author. W 
the further increase in discharge, the flow becomes completely turbulent, a 
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yin the roll waves are predominant in forms of bore if the channel rough- 
3S is approximately less than 0.01 in Manning’s roughness coefficient, 
ereas for rough channels in which the coefficient is larger than 0.01, the 
I waves can not be observed. When the discharge is still increased, the roll 
ves ultimately appear even in such rough channels. This fact is evidently 
lal to the results obtained by Thomas, 10 Dressler]! and others. 

In Fig. 1, our experimental data for roll waves are plotted, adding supple- 
ntal data of the author and Binnie. It is then concluded that the foregoing 
atment for the formation of roll waves is verified by the experimental ob- 
vations. Evidently, the law of resistance and the momentum correction 
tor play a very important role in the criterion for formation of roll waves. 
. 2 indicates the criterion in terms of velocity formula and momentum 
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‘IG. 2.-ILLUSTRATION OF CRITERION FOR FORMATION OF ROLL WAVES IN 
TERMS OF LAW OF RESISTANCE AND MOMENTUM CORRECTION FACTOR. 


rection factor, and it also illustrates that roll waves cannot be observed in 
transition region from laminar to turbulent flows, as (1 + b)/a is very 


all in this region. : é 
In the same manner, the criterion for roll wave formation in triangular 


nnels will be obtained as follows. 


N2, = 4/sin $ | (222)? - 4(- n(222) - a(e- of ee (3f) 


! . Fig. 3 indicates the 
yhich ¢ is a half of the angle between two side walls 

1e criterion for triangular channels, and it is obvious that the roll-wave 
mation in triangular channels cannot be compared with that in rectangular 
nnels. 

[he behavior of wavevelocity of roll waves when the wave pattern pesaures 
mately stable will be described next. By the use of the expressions for 
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progressive discharge rate and for the mean velocity of flow at the cont: 
section, the wave velocity is determined. Introducing the dimensionless wa 
velocity by Vw' = V. w/t do Cos 9)1/3 for recede channels, in which qg 
the discharge per unit width for uniform flow, Vy,' 


Vw' ={a- 1+ Ya(a- 1) + TREY es PAA Ores (3 


As the dimensionless progressive discharge rate K', in which 


Ki = K & cos 8 


Ww 


1 I Va(a@ - 1) + eu a - Va(a- 1) + ey ( 
x2 oe 


a - (1/N2) a - (1/N%) 


is a function of a and Np. so the dimensionless wave velocity is also é 
pressed in terms of the momentum correction factor and Froude numbe 
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FIG. 3.—CRITERION FOR FORMATION OF ROLL WAVES IN TRIANGULAR 
CHANNELS. 


Fig. 4 shows the behaviors of Vy' for both laminar and turbulent flows. In 1 
same figure, our experimental data are plotted with data of the author a 
Binnie, and it is found that the experimental verification, to the theoretic 
approach, can be made. 

By the similar procedures, the behaviors of wave length, wave period, 2 
wave height as the hydraulic characteristics are determined. Introducing 1 
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nd dimensionless parameters of wave length, wave period and wave 
ight, 


eer easing te GOs ey Va jens eee ake it (34a) 
EAL Gen aarp aeaiedeeaet pg pian (34b) 
d Dj, = Dp (g cos Bevo), Ore peor abies (34c) 


gs. 5-7 are graphically obtained by the numerical calculation, using the 
regoing relationships. The actual wave height at the frontal side is always 
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FIG. 4.—BEHAVIORS OF DIMENSIONLESS WAVE VELOCITY. 


s than the calculated one, because the calculated wave profiles have sharp 
ve crests due to the discontinuous shock front, whereas the actual wave 
st is round probably because of the capillary effect and the existence of 
ticities at the frontal side. 

Other hydraulic characteristics in flow, due to the formation of roll waves, 
. calculated in the same way. However, these are not concerned herein, and 
ees of hydraulic characteristics in flow will be seen in the writers’ 
er. *~ 

Based on the results derived and the knowledge obtained by the writers 
ough the research works, the following 16 items are presented: 

1. The author states that the roll wave is the result of interaction between 
surface tension and the gravity force in a slightly disturbed flow. How- 
2 


‘ 
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ever, the hydraulic characteristics of roll waves in laminar flow can be thec 
retically analyzed as presented in the foregoing, without introducing the effer 
of surface tensions, and the results derived theoretically are satisfactorily 1 
agreement with the experimental ones obtained by the writers and Binni 
From this fact it may be concluded that the effect of capillary forces on tt 
formation and the characteristics of roll waves is quite little or rather negli 
gible practically, except in cases of extremely steep slopes and of very lo 
Reynolds number. 

2. The author concludes that slug flows result from instability whic 
causes the transition from supercritical laminar to turbulent flow. Slug flow: 
defined by the author, seem to be corresponding to incomplete roll waves : 
low Reynolds number in turbulent flow just changed from laminar to con 
pletely turbulent through transition in its flow regime. The fact that roll wave 
and slug flows are not observed at the transition from laminar to turbule 
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FIG. 5.—RELATION BETWEEN DIMENSIONLESS WAVE LENGTH AND FROUDE 
NUMBER. 


flow is caused by the different nature of flow and thus the different laws 
resistance as described in the foregoing. Slug flows, therefore, should be e: 
plained as kinds of roll waves which result from instabilities of flow itself, 
low Reynolds number in turbulent flow as well as roll waves in laminar flow 

3. The author’s Eq. 11 is questionable, whether it represents the criteri 
for formation of roll waves and slug flows or not. According to Eq. 11, rc 
waves are formed in laminar flow, when Ny <2. However, the author’s Fi 
12 indicates that roll waves were ahapevell : even when Np > 2, and the write 
also obtained the same results as shown in the figures uateated in the for 
going. Consequently, Np = 2 is not a criterion for formation of roll waves 
laminar flow. 

4. The author applied the criterion for the instability of turbulent flow a 
the formation of turbulent roll waves, : 
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hich was derived by the use of the Chezy formula, to laminar flow. Since, 
wever, Eq. 35 was derived by putting C = constant, the author’s Eq. 21 
ould not be used to obtain Eq. 17a from Eq. 35. 

5. In regard to the effect of surface tensions, the author suggests that the 
locity of roll waves is expressed by 


Bape et Cate orci alee dawndany (36) 
1is relation is based on the idea that the wave length decreases and con- 
rsely the wave height increases with time, and ultimately the condition of 
inimum celerity C,,;, is attained. However, it is evident from Figs. 20 and 
| that the wave velocity, the wave length, and the wave height actually in- 
‘ease until the terminal values are reached with the lapse of time. This fact 
inconsistent with the author’s suggestion. 
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FIG. 6.—RELATION BETWEEN DIMENSIONLESS WAVE PERIOD AND FROUDE 
NUMBER. 


6. For laminar flow, the author’s Eq. 19 verifies that the plot of Np/ g1/ é 
ainst Rp shown in Fig. 14 is correct as illustrated again in Fig. 18. In the 
me manner, since Eq. 19 is written as 


N/VS = const. (g, v, o/p) - (Nw/$1/3)3/10 Fuipca: Cae (37) 


e author’s Eq. 15 is suggested to plot N,/S1/2 against Ny/s!/3. . 

7. Graphical expressions for Np against Rp in Fig. 12 and for Np against 
, in Fig. 13 are the indication of the laws of resistance in each regime of 
w, and the hydraulic characteristics of roll waves and slug flows are not 
wn. Therefore, the author’s conclusion is not derived from these figures, 
t the surface tensions and viscous effects predominate in the phenomenon 
roll wave formation, and at higher Reynolds numbers, corresponding to the 
ditions which led to slug flows, viscous and surface tension effects de- 
~ased and, finally, became negligible. 

According to the results of analysis by the writers, the viscous effect on 
characteristics of roll waves in laminar flow is introduced indirectly as 
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Froude number from the relation of Eq. 19 because all the dimensionle: 
quantities of roll wave characteristics are functions of the Froude number, : 
shown previously. On the other hand, characteristics of roll waves (slug flow 
in turbulent flow are not influenced by the viscosity as far as Manning 
formula is used and its roughness coefficient is assumed independent on tl 
viscosity. Actually, however, the viscous effect may exist more or less if tl 
channel bed condition is a hydraulically smooth boundary. 


8. If Np and For are put in a form of 


it is found after simple computation using the author’s Eqs. 19 and 37 that fl 
value of 6 should be 1/4 for laminar flow instead of 0.218 in the author’s E 
32. 
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_ FIG. 7.—RELATION BETWEEN DIMENSIONLESS WAVE HEIGHT AND FROUDE 
NUMBER. 


9. In obtaining the minimum slope 0.029 ~ apc for the formation of r 
waves, the criterion to laminar flow 
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was applied. This criterion is corresponding to 
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wever, Fig. 17 indicates that roll waves develop even when NF < 2 and it is 
Gent that the latter is inconsistent with the author’s idea for the roll wave 
‘mation. 

10. In deriving the minimum slope 0.01 for the formation of slug flows, 
. 39 was also applied. At the Reynolds number 1,200, which the author used 
a lower limit of Reynolds numbers, the flow regime will be the transition 
m laminar to turbulent flow. Therefore, it is doubtful that Eq. 39, which 
s derived by using the law of resistance to laminar flow, can be applied to 
‘formation of slug flows. 

11. The writers are much interested in Fig. 19 which shows the location of 
tiation of roll waves and slug flows. This figure indicates the fact that roll 
ves and slug flows cannot be formed in the transition from laminar to tur- 
ent flow, otherwise very long distances are needed for their development. 
is fact was already pointed out by the writers in their Fig. 1. 
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}. 8.-PLOTS OF THE WRITERS’ AND BINNIE’S DATA OF ROLL WAVE VELOCITY 
IN FORM PROPOSED BY THE AUTHOR. 


12. The author suggests that the wave velocity of roll waves Vw is equal to 
-Cmin» in which C,pjn = 0-763 ft per sec. Although the value of Vy - V ob- 
1ed from the experimental data is of the same order in magnitude as Cin, 
suggestion to the wave velocity is too intuitive and the effect of capillary 
siven too much evaluation. pals 

13. The writer’s and Binnie’s data of the wave velocity in laminar flow are 
sented in Fig. 8 in the same way as plotted in Fig. 24. Itis dound from the 
ire that the author’s data are larger than the writer’s and Binnie s, except 
Reynolds numbers of 280 to 420. This disagreement in addition to that in 
. 4 cannot be understood by the writers. 

14. The ratio of the wave velocity Vy to the surface velocity Vg, in laminar 
7 is written, from Eq. 8, as 


a 
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This relation shows that the ratio Vy/Vg is a function not only of Rp but al: 
of the slope S. Therefore, the slope should be introduced as a parameter 
Fig. 25. It is not difficult to find from Eq. 42 that the value of Vy/V is alway 
greater than unity and increases with decrease in Rp if the slope is constan 
This trend agrees with that shown in Fig. 25. 

15. In turbulent flow, the ratio of the wave velocity (slug flows) to the su: 
face velocity is expressed, from Eq. 6 if assuming a@ = 1, as 


%-(E) (+8) « Ge) (rks) 


Although it seems possible from Eq. 43 that the value of Vw/V, becomes le: 
than unity, which depends on the value of V/Vg, the results of his observatio 
which indicates that the surface velocity exceeds the wave velocity of sh 
flows, can not be understood. If so, when the flow with waves is transferred 
steady-state by moving the coordinate system at a constant wave velocity, t 
water surface in the virtual flow of water running upstream has to have 
velocity in the downstream direction. This phenomenon is possible only in t 
vicinity of the wave crest which results from the horizontal vortex due 
hydraulic jump. 

By only this possible behavior of flow, the experimental results of sl 
flows cannot be interpretated analytically. 

16. It may be rather better than the plot in Fig. 26, that the ratio of ft 
wave height to the wave length D}/L is plotted against Froude number with 
parameter of the slope because both Dp g cos 6/ Veg" and L tan 6g cos a/V, 
are functions of Np as presented previously. 


SPILLWAY DESIGN FOR PACIFIC NORTHWEST PROJECTS2 
Discussion by B. Michel and A.R. Gagnon 


B. MICHEL! and A.R. GAGNON,?2 A.M. ASCE.—This paper has brought up 
2 of the oldest, but still one of the most discussed, problems of hydraulic 
yineering: the design of spillway crests. 

The author should be commended for having presented his interesting views 
this problem and for the attempt to standardize the design of spillways. 

In this last regard the writers agree fully with Mr. Webster that semicir- 
ar pier noses induce small flow contractions and that a 1 on 12 taper on the 
ynstream end of gate slots is generally one of the best means to lower cavi- 
ion risks at the slots. 

However, the writers believe that a second thought should be given before 
eralizing the process of designing spillway cress for 75% of the nominal 
d. 

An outstanding contribution in that field has been made recently by R. 
noine? and the writers believe that he has given a very general basis for 
llway crest design. Mr. Lemoine proposes to adopt a cavitation factor for 
llway crests in the same manner as for the design of turbines: 

This factor is 


vhich p,, is the water vapor pressure at the site, p represents the absolute 
imum pressure on the spillway crest, v2/2 g is the velocity head at the 
at of minimum pressure, and w denotes the specific weight of water. 

This factor takes into account, not only the actual permanent negative 
ssure in the dangerousarea, but also thedanger of having the velocity head 
isformed into an additional negative pressure in low pressure nuclei ed- 
3, formed behind asperities of the spillway surface. Thus, the factor should 
er all conditions of cavitation on spillway profiles. | 

[he design. of any spillway can then be made with the following considera- 


r= 0 : certain cavitation; 


Dee — ee ee 
August, 1959, by Marvin J. Webster. 

Prof. of Hydr. Engrg., Laval Univ. 

‘Lecturer in Fluid Mech., Laval Univ., Quebec City. 
“Cavitation on Weir Crests,” by R. Lemoine. Procee 
R, The Hague, 1955. 5 
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o = 1 : cavitation probably impossible, even if all the velocity head w 
transformed into additional negative pressure; and 

o = 0.5 : acceptable risk of cavitation if this working condition occu 
rarely and if the spillway surface is very smooth. 


For a standard ogee spillway it is easy to know the minimum pressure 
the crest,4,5,6 the velocity head, and to design the spillway with the help 
the cavitation factor. It is indeed interesting to note that the foregoing refe 
ences show that, for an overdesign of 75%, the cavitation factor, after a “c1 
tical design head,” is invariably zero while the minimum pressure becom 
equal to the vapor pressure. It might be also noticed that the classical fol 
is not generally the one that will give the highest discharge coefficient fo1 
certain value of the cavitation factor. Model tests must be made in each cas 

The model results given in Mr. Webster’s paper, for the Chief Jose 
Spillway, show a minimum pressure of -23 ft of water with no piers and -9 
with piers. These results are somewhat higher than those obtained on simil 
spillway profiles by many investigators. This could probably be explained 
the fact that very small differences in the geometry of each project have 
important effect on pressure values, and also by the very high difficulty 1 
volved in measuring the peak of rapidly varied pressure distribution on 
scale model. The writers do not believe that the results for this spillway a 
representative of all those of the same type. 

Taking the last value for the minimum pressure of -9 ft with an operati 
head of 54 ft it is found that the cavitation factor would then be 0.38, which 
somewhat lower than could be permissible with Mr. Lemoine’s theory. 

The author was fully aware of some cavitation risks and states that due 
the low predicted frequency of the maximum discharge, the cost saved in i 
itial construction more than offsets the cost of infrequent repairs resulti 
from damage due to negative pressure. 

The writers agree with him and this could easily be taken into account 
choosing, in each case, anaccepted value of the cavitation factor. The write 
strongly believe that this would be more logical than to accept a percenta 
head design irrespective of the value of the head itself, the actual shape 2 
pressures on the crest. 


> 


4 “Untersuchugen an Uberfallen,” by O. Dillmann. Mitteilungen des hydrate 
Institus, Munich #7, 1933. 

5 “Engineering Hydraulics,” Hunter Rouse (Editor), Chapter VIII. 

6 “Model Research on Spillway Crest,” by Rouse and Reid. Civ. Enero vol 


No. 1, 1935. 


FLOOD CONTROL ASPECTS OF CAUCA VALLEY DEVELOPMENT2 
Discussion by Julio Escobar- Fernandez 


JULIO ESCOBAR-FERNANDEZ.!—Three aspects from the paper by 
ssrs. Kirpich and Ospina could be considered. First, the significance and 
ortance, in the development of the so-called underdeveloped countries, of 
establishment of regional autonomous corporations of development; sec- 
, the technical difficulties that have to be faced by engineers interested in 
er problems in such countries, where hydrological data is invariably 
rce, or of doubtful value and inconsistent; and third, the technical features 
big hydroelectric project, fleetingly mentioned in the paper, which could 
developed by diverting the waters of the Cauca River to the Pacific coast, 
ystem that could be classified as one of the important energetic possibili- 
, not only of Latin America but of the world, with the particular advantage 
ts proximity to the sea, a feature that makes it especially attractive for 
establishment of heavy electrochemical industries. 
\lthough the three aspects might be of general interest, due to limitations 
pace, only the first one will be discussed, at present. Some of the general 
ciples governing the policies and the organization of regional corporations 
sted to local development are set forth, and some emphasis is placed on 
importance of such bodies as helpful instruments in the execution of wider 
lopment programs, usually national in scope. The timeliness of the sub- 
seems especially relevant today (1960), when better management for 
1omic improvement constitutes a daily concern for the governments in the 
; developed countries. 
‘olombia is a country with rich natural resources that can be developed to 
ite an economically solid nation, thereby raising the standard of living of 
eople. To reach these two objectives, in time and without waste of effort, 
; considered that a plan for the development and effective use of those re- 
rces is indispensable, as is the implementation of the plan. In Colombia, 
Central Government has already launched such a program and the recent- 
rganized Planning Unit, directly responsible to the President of the Re- 
ic, is a fundamental step towards the coordinated execution of the plans 
have been already formulated or which are under study. Asa starting 
t, a better orientation of government investments is considered one of the 
ary policies of development, because such investments as in transporta- 
. energy production, irrigation and public utilities are or should be the 
dation on which a solid economy must rest. 
m the other hand, the expansion of public enterprise, which constitutes the 
r fundamental basis on which the progress of nations depends, must be 
weet See ee 


‘September, 1959, by Phillip Z. Kirpich and Carlos S. Ospina. 
‘Escobar, Venegas y Rodriguez, Ingenieros Ltda., Bogota, Colombia. 
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stimulated not only by the facilities provided by those public investments, k 
by the creation, through government measures, of a minimum of general co 
ditions favoring private investment. 

As one of its immediate tasks, the new Planning Unit will therefore have 
establisha balance between all resources available—economical, physical a 
human—and the many needs of the country. From this overall picture, a mo 
or less rigid schedule of priorities must be determined, in order to establi 
how and where such resources can be spent most usefully and effectively. 

Once a general plan for the whole country has been defined, the next st 
should be to project it on a sectorial basis, to use the economist’s jargon, 
on a regional basis, if the characteristics of the country so allow, as ist 
case in Colombia. 

A program of integrated regional development, such as the one undertak 
in Colombia in the valley of the Cauca River has, among others, the followi 
advantages: 


1. It permits the establishment of a kind of experimental field where se 
eral techniques of development can be tested. In this way, the selected zo 
becomes a demonstration area, where direct experience is the best school; 

2. financial problems may be of a reduced magnitude in comparison witl 
nation-wide program of development. 

3. administrative problems are also reduced; and, 

4. in developing a particular area, the economic benefits spread thoroug 
out the country. 


When in 1954, David Lilienthal, mentioned by the authors of the paper und 
discussion, recommended to the Government of Colombia the creation of 
autonomous regional authority,2 as a first stage in its wider program of d 
velopment, he based his recommendations on several general principles tl 
should guide this type of body, principles that can be resumed as follows: 

The initiative for a general program must originate in the region itself. 
is therefore indispensable that people in the area be willing to cooperate 
the plans, and that enthusiasm for innovation and improvement be general. 

Second, the approval and cooperation of the Central Government is al 
necessary. The Government should be disposed to frame the program w 
some autonomy, through the creation of an agency or public corporation ¢ 
voted to that purpose. This agency should be provided with sufficient admi 
istrative authority, and with well defined functions and aims. To better defi 
the status of that corporation, it may be useful to quote the words of Presid 
Roosevelt when he asked Congress to create the TVA in 1933 (quoted by N 
Lilienthal,: The authority should be, he said “a corporation clothed with 
power of government but possessed of the flexibility and initiative of a priv: 
enterprise.” 

Third, the financial assistance of the Government or of any internatioi 
organization, such as the World Bank. 

Fourth, the creation of a local economic development corporation to p1 
mote and stimulate private enterprise in its own development programs. 

Finally, it is advisable that the selected region for an experiment of tl 
sort, especially if it is the first to be undertaken in the country, should be 


2 “Recommendations on the Establishment of Regional Development Authorities 
the Republic of Colombia,” by David E. Lilientahl, June 25, 1954. 
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re the benefits obtained through the techniques of planning, are more or 
$s immediate, in order to overcome the natural tendency of people to be 
ptical or even hostile to programs that require a long time to produce 
nite results. 

All of these conditions were satisfied by the Colombian region of the Cauca 
ley, and the CVC, the autonomous corporation created in 1955, is beginning 
show results in a way that may be considered a real success. 


FRICTION FACTORS IN CORRUGATED METAL PIPE? 


Discussion by A.R. Chamberlain and Nicholas Bilonok 


A.R. CHAMBERLAIN, ! A.M. ASCE.—The authors are to be commended for 
ir very careful work in obtaining design information on friction factors for 
ze-diameter corrugated pipes. These data and results presented by the 
n0rs will undoubtedly be used by design engineers for several decades, in 
cases where pipes of thenature tested are used. This paper represents the 
ult of a tremendous amount of effort on the part of Army Corps personnel. 
The writer wishes to extend the data presented by the authors on the basis 
slear water, full pipe flow data obtained during 1953-56 in the Hydraulics 
oratory, Colorado State University, Fort Collins, Colo. These data were 
uined incidental to certain basic research studies on nominally 12-in. di- 
ster smooth, standard corrugated and helical corrugated pipes. 
Zach of the types of corrugated metal pipe sections tested were 100-ft long, 
on a zero slope. Both corrugated pipes had a clear inside diameter of 
| + 0.03 in., as determined from ten trials on each pipe. 
[The water was recirculated in a closed system. The discharge was meas- 
d by means of a ’sharp-edged orifice which had been calibrated volumetri- 
y while inplace. Discharges were determined, for eachdischarge-control- 
fe setting, from the orifice calibration curve using the arithmetic average 
our separate readings of the differential head across the orifice. 
lydraulic gradient was determined by means of water piezometers located 
0-ft intervals along the pipeline. The piezometer openings into the pipeline 
e 3/64-in. indiameter. For each discharge, each piezometer was read 
‘ times and the arithmetic average used in hydraulic gradient computa- 
s. Upstream end (for 20 pipe diameters) and downstream end (for 10 pipe 
neters) piezometers were generally ignored, though recorded, in hydrau- 
sradient studies. The piezometers were located on the crest of the cor- 
itions (viewer inside the pipe) after preliminary investigations of the ef- 
of piezometer location along the corrugation, comparable to those of the 
or’s. 
‘he procedure was to first read the orifice manometer, then the piezome- 
sank and next the orifice manometer again, until each had been read four 
‘s during a period of about lhr. Next the discharge valve setting was 
ged, 15 to 20 min. allowed for the flow to stabilize and the next set of 
to be obtained. / 
‘ig. C1 shows the two boundaries tested and reported in this discussion. 
standard corrugated pipe was zinc coated close-riveted construction. 


September, 1959, by M. J. Webster and L. R. Metcalf. : 
ting Dean, Coll. of Engrg. and Chf., Engrg. Research, Colorado State Univ., Fort 
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The corrugations had a pitch of 0.222 ft along the axis of the pipe. The amp 
tude of the corrugations from the mean was 0.0208 ft. The crests and trous 
were circular arcs with an included angle of approximately 53.73 deg anc 
radius of 0.0573 ft. Straight tangent sections joined the circular arcs. 

The helical corrugated pipe was zinc coated. It had a continuous lock se 
‘ joint. The corrugations had a pitch or wave length of 0.167 ft normal to” 
axis of the corrugation. The amplitude of the helical corrugations was 0.0: 
ft, measured from the mean elevation. To complete the corrugations, circu 
arcs witha radius of 0.0370 ft were joined by straight tangent sections. 14 


) 
) 


12-in. Standard Corrugated Pipe 
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\2-in Helcor Pipe 


Note: 
Alllength dimensions are in inches. 
Section views are normal to corrugation. 


FIG. C1.—DETAILS OF BOUNDARY FORMS, 


helix angle, measured between a line drawn along the outer extremity of 
pipe parallel to the pipe axis and a line parallel to a corrugation, was 66 de 

Garde’s data were obtained with clear water full pipe flow after the pi 
lines had operated several hundred hours, transporting 0.2mm mean diame 
sand at varying concentrations. Chamberlain’s data, designated by a cross. 
the Helcor pipe only, were taken just prior to Garde’s. The remaining d: 
all taken by Chamberlain, were recorded with clear water full pipe flows 
new pipe. 

Derived quantities Re, f, and V were either computed directly or read f1 
graphs prepared for the purpose. In either case, the pipe diameter emplo 
in the analysis was 1.00 ft, slightly less than the inside clear diameter. 1 
magnitude of the diameter was used because (a) it was very convenient, 
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re important (b) it was felt that most engineers applying the data would 
omatically use 1.00 ft in all their design computations. 

Tables 1 and 2 summarize the full pipe flow data obtained, admittedly not 
*y extensive, but hopefully of some use to the profession. These data are 
sented in Fig. C2, for both the standard corrugated and helical corrugated, 
a manner comparable to the 
hor’s Fig. 10. They are also 
en in Fig. C3 in the same way 10.0 
the author’s Fig. 11. ao 
Fig. C3 indicates that a rec- 
mendation for new 1-ft diam- 
r helical corrugated pipe, of 
0.040, is appropriate. This 
ne figure shows that f is nota 
ction of Reynold’s number 
r the range tested. This may 
partially explained by the fact 
t the corrugations are not 


Slope of Hydraulic Grade Line in Per Cent 
9° 
@ 
te) 


‘mal to the direction of flow in 1.0 
pipe. 
The writer is unable to say 060 
, Garde’s data, previously un- o.soL 
lished, have such a wide scat- 0.40 
This is particularly true of 
helical corrugated pipe data, 0.30 
ce the pipeline was not in use 0.20 
ing the interval between : 
mberlain’s and Garde’s used- © Garde 
e tests. : 
Perhaps the most important 0.10 sarees BS 
at to note is that the dashed < Chambarialn 
tion at the smaller diameter 
of the author’s Fig. 17 can be 0.05 


“0.5 O08! 2 3 45678910 


anded by plotting in for stand- Discharge in CFS. 


corrugated pipe the tentative- 
ecommended point f=0.12,at = rq. c2.—DISCHARGE VERSUS HYDRAULIC 
1.0 ft. Alsoto be plotted on = GRADIENT FOR 1-FT. DIAMETER HELCOR 
. 17 for corrugated pipe is the AND CORRUGATED PIPES. 

puted (from f = 0.120) 
0.0255 at a pipe diameter of 1 ft. The data are insufficient to determine if 
a function of Reynold’s number. 

‘ig. C4 graphically describes the variation of the total differential head 
veen crest and trough piezometers for the standard corrugated pipe tests 
. function of Reynold’s number. Perhaps more important, however, is the 
, also in Fig. C4, of the ratio of the hydraulic gradients as determined hy 
meters at the crest(s) and at the trough(s) of the corrugations, as viewed 
n inside the pipeline. This plot substantiates that the author’s and the 
er’s method of determining hydraulic gradient from a set of piezometers, 
ted consistently at the same relative location on the corrugations along 
pipeline, is sufficiently accurate for engineering purposes over the range 
eynold’s number incorporated in the tests. 
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TABLE 1.—SUMMARY OF TEST RESULTS FULL PIPE FLOW 
1-FT-DIAMETER CORRUGATED PIPE# 


V, 
in feet 
per second 


(1) 


Water 
Temperature 
in C° 


Q, 
in cubic feet 
per second 


Re 
x 1073 


: 590 0. 
6.654 583 0. 
7.26% 655 0. 
7.264 660 0. 
g 
2.38b 0. 
4.45b 0. 
5.30b 0 
5.98b 0. 
7.23 0. 
0. 


2 Data from R.J. Garde, “Sediment Transport Through Pipes,” Thesis in partial 
fillment of the requirements for the Master of Science Degree, Colorado State Univ 
sity, Fort Collins, Colorado, December 1956. 

Data from A.R. Chamberlain, “Effect of Boundary Form on Fine Sand Transpor 
Twelve-Inch Pipes,” Dissertation in partial fulfillment of the requirements for 


Doctor of Philosophy Degree, Colorado State University, Fort Collins, Colorado, J 
1955. 


[JI FT. DIAM. STO. 
CORRUGATED 


qovaclay Number "Re" 


Darcy- Weisbach "f" 


FIG. C3.—DARCY-WEISBACH COEFFICIENT VERSUS REYNOLD’S NUMBER FOE 
1 FT, DIAMETER HELCOR AND STANDARD CORRUGATED PIPES. 
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TABLE 2.-SUMMARY OF TEST RESULTS FULL PIPE FLOW 
1-FT-DIAMETER HELICAL CORRUGATED PIPE 


Vv, Q, Water 
n feet in cubic feet Ss, Temperature Re 
r second per second in % in °C x 10-3 f 

(1) (2) (3) (4) (5) (6) 
1.592 1.25 0.25 16.5 135 0.0640 
2.628 2.06 0.50 15.8 218 0.0478 
2.76 27 0.59 16.7 235 0.0500 
3.314 2.60 0.95 17.0 284 0.0559 
3.582 2.81 0.96 16.2 302 0.0483 
4,174 3.27 1.27 16.7 356 0.0470 
5.732 4.50 2.20 17.0 493 0.0435 
7.502 5.89 3.56 14.6 608 0.0408 
7.978 6.25 4.02 15.9 666 0.0408 
8.35% 6.56 4,22 16.3 705 0.0390 
0.98 0:75 0.062 15.2 82 0.0413 
1.055 0.80 0.068 15.4 89 0.0396 
1.10 0.84 0.074 15.0 91 0.0394 
1,46) 0.0368 
1.63 0.0401 
1.772 0.0390 
1.88? 0.0375 
2.08» 0.0446 
2.165 0.0397 
2.22 0.0412 
2.45 0.0400 
3.43 0.0384 
4,10 0.0360 
4.29 0.0385 
4.45¢ 0.0422 
4,520 0.0404 
4,80b 0.0419 
5.17 0.0391 
5.41b 0.0382 
5.45) 0.0412 
5.80) 0.0383 
6.00° 0.0382 
6.23 0.0397 
6.61 0.0391 
6.80 0.0390 
7.23¢ 0.0418 
7 94b 0.0397 
7.39¢ 0.0350 
7.40b 0.0376 
795b 0.0433 
s.10b 0.0422 
3.210 0.0411 
8.40¢ 0.0406 
ane 0.0392 
B.46 0.0396 

gab 0.0409 
cP 0.0396 
9.51) 0.0392 
9.72 0.0456 
9.80¢ 0.0438 


1 Data from R.J. Garde, “Sediment Transport Through Pipes,” Thesis in partial fulfillment of the 
uirements for the Master of Science Degree, Colorado State University, Fort Collins, Colorado, 
ee Chamberlain, “Effect of Boundary Form on Fine Sand Transport in Twelve-Inch 
es,” Dissertation in partial fulfillment of the requirements for the Doctor of Philosophy Degree, 
gee State University, Fort Collins, Colorado, June 1955.. 

‘Data obtained by Chamberlain after pipeline had been used to transpo 
red hours. 

Arithmetic average, all Chamberlain data, f = 0.0400. 


rt 0.2 mm sand for several 


72 March, 1960 HY 


Acknowledgements.—The basic research studies were supported by Arm« 
Drainage and Metal Products, Inc., Tau Beta Pi Association (through a fellov 
ship), and Colorado State University. 


H|n 1.0 


0.9 


° 
@ 


9° 
a) 


° 
b 


Diff. Head, Trough to Crest (ft) 


9° 
nS) 


FIG. C4.—RATIO OF HYDRAULIC GRADIENT ON CORRUGATION HIGHS TO THAT C 
CORRUGATION LOWS, AND DIFFERENCE IN HEAD BETWEEN HIGHS AND LOV 
VERSUS REYNOLDS NUMBER FOR 12-IN. CORRUGATED PIPE. 


NICHOLAS BILONOK.2—The most important contribution by the authors 
the presentation of data from certain carefully conducted experiments on fi 
size, large corrugated metal pipes, to determine friction factor “f” and abs 
lute roughness “n.” The accuracy depends on a correct evaluation of “n,” t 
well known coefficient of absolute roughness in the Manning formula, whi 
depends primarily on the surface condition of the perimeter of a pipe or ope 
channel, 

Because of the complex nature of the friction factor “f” various attemp 
have been made to group all the roughness dimensions together into a sing 
coefficient of absolute roughness “n” for simplicity. , 3 

The writer is particularly interested in this paper as he had independent 
investigated the relationship between friction factor “f” and absolute roug 
ness “n,” for many years. | 


2 Civ. Engr., Engrg. Div., Corps of Engrs., U.S. Army, Philadelphia, Pa. 
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The values of Q, S, Hg, L and D have been measured in experiments and the 
iction factor “f” has been found to be given by 


ne equation for the roughness coefficient “n” for the Manning formula is 


1.486 R2/3 g1/2 
n = Sra te ee i he Cae ee (3) 


j. 3 may be written from the friction factor “f” by the relationship 


1. 4 can be written as 


_ 116.57 n2 
R1/3 


‘substituting R = D/4 


‘om a review experiments data presented in the paper the writer has com- 
led Table 1 to relate the friction factor “f” to Manning’s roughness coeffi- 
ants “n” ranging from 0.010 to 0.025. 

The values of “f” computed from Eq. 6 are in agreement with those ob- 
ned by experiments for pipes 3 ft, 5 ft, and 7 ft in diameter, and also with 
ble V in Reference 3. ; 

Another quantity of importance that is the Reynold’s number: 


The exact form and position of “f” versus Reynold’s number for a given 
Ye diameter and roughness is not clear (Fig. 11), and close correlation re- 
ires more study. 

The practical value of the Reynolds number is that it indicates the degree 
turbulence in flowing water. For a given size of circular pipe the velocity 
the major variable and the Reynolds number increases as the velocity of 
w increases. 


: 
7 
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The author’s Tables A and B indicate that the Reynolds number varies 
velocity of flow only. The water temperature changes have no significance. 

Investigations of friction losses for Part-Full flow, in Reference 3, ha 
found that “For uniform subcritical flow, the Manning coefficient did not inc 


TABLE 1.—VALUES OF f x 102 


0.011 0.012 0.013 0.015 0.018 0.021 0.024 


DBmamuouwndod 
RPwWOON IN Oo 


cate a systematic variation with the Reynolds number or pipe size within tl 
accuracy of experimental data.” 

It can be recognized that several factors might have affected the accurac 
of the test results and the computations. 


REVISED COMPUTATION OF A VELOCITY HEAD WEIGHTED VALUE® 


Discussion by Steponas Kolupaila, Israel H. Steinberg 
and William C. Peterson 


STEPONAS KOLUPAILA.!—Presentation of this paper proves that the 
yblem of the correction factor, known as a Coriolis’ coefficient, is still 
ive, and still not entirely clear. The purpose of this factor is to correct 
‘velocity head in the Bernoulli equation for different velocities across the 
tion. 

The basic expression given under 10, p. 72 is correct, except for an ob- 
us typographic omission: AQ is to be understood instead of Q in the nu- 
rator. This formulacan give an accurate solution when a stream is divided 
9a large number of small areas with measured velocities. This work is 
re reasonably performed by a graphical integration with a planimeter, 
minating those “astronomical” numbers. Application of this formula for 
reral parts of a cross section with assumed coefficients is a half measure. 
= value obtained by the authors (1.19) seems to be too low for anopen chan- 
with overflooded banks.2 The river gagings by two point or one point 
thods, although more than sufficient for the routine discharge measure- 
nt, are positively inadequate for the determination of the Coriolis coeffi- 
nt. Two integrations are necessary in an open channel, along the verticals, 
| across the channel, and more points must be taken on every vertical. 


ISRAEL H. STEINBERG,? F. ASCE.—The sample computations presented 
he article, and which are patterned after a procedure given in the Engi- 
ring Manual of the Corps of Engineers, still contain some of the objection- 
e, relatively large, figures required to obtain a relatively small value of 
ocity head. For this reason the writer has adopted a number of alternate 
cedures in order to further reduce the necessity of dealing with large and 
rieldly values. One of the methods utilizes a weighted area curve for the 
tion. The equation is: 


1 
Ag = Ae) eee ayer (1) 
a! (ania) 2| 1/2 
Cc c Ag 


Meee ie OS 


. September, 1959, by J.M. Lara and K.B. Schroeder. 

L Prof. of Civ. Engrg., Univ. of Notre Dame, Notre Dame, Ind. 

: “Open-Channel Hydraulics,” by V.T. Chow. McGraw-Hill, 1959, p. 28. , 

} Asst. Chf., Planning and Reports Branch, U.S. Army Engr. Dist., San Francisco, 


ps of Engrs. 


a 
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in which Aw is the weighted area, Ac and Ao are respectively, the channel a 
overbank areas, and Cc and Cy are, respectively, the channel and overba 
conveyance factors of the form: 


1.486 AR2/3 kK 
Peis oh eels why 


eoeeaeveses n8 6s oe 


in which C is the total conveyance factor for the section (=Ce +Co), and L 
the reach length. 

Eq. 1 was developed for flow lineproblems involving separation of the se 
tion into two parts; namely, channel and overbank. The equationcan be read 
expanded to include any desired number of components. Taking the reciproc 
of Eq. 1 and squaring gives: 


(a) = (ae) UCB) (2G) 
POuGueChesc 


Eq. 3 applies, also, for any number of components and can be written in t 
more generalized form of 


(z) = a)". 2) CP 
alas eerie | 


If L is constant throughout, the K-values can be substituted for correspondi 
C-values. 

Eq. 4 can be further modified in more convenient forms for computi 
either (a) the weighted velocity head, hy, or (b) the alpha coefficient, a, 
terms of total section area. To obtain hy, both sides of Eq. 4 are multipli 
by (Q/8.03)2, resulting in 


( Q Ns (es Q iE [(c2 3/2 Q 2 
Shes oN AOS Rega Moe. 8.03 Aj * (2) 8.03 7 | 


Application of Eq. 5 in computing the weighted velocity head is illustrat 
in Table 1 using the author’s values for section 2. The values in Cols. 7 and 
were obtained with simple settings of the slide rule. 
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TABLE 1 


66,400 0.255 PRE: 1.511 0.149 
67,000 0.258 2.295 0.352 

5,240 
26,300 


ff the value of a is desired, Eq. 4 is multiplied by A2, the total sectional 
a, to give: 


Boe. +1: 
‘Zea 


|. 


yomputations are the same as for the 
mple above except for Cols. 8 and 9. 
[Table 2 only Cols. 7, 8 and 9 are 
vn for brevity. 

‘or the assumed value of Q = 10,800, 


+ 


+ 


10800 
hy = 1.424 (a3 X 2226 


)* = 0.521 


VILLIAM C. PETERSON, #4 M. ASCE.—Computation of an approximation of 
velocity head correction factor, a, by the fraction 


2 
K3/a2 + K3/a2 peice K3/A2 
(2K)3/(ZA)2 


ssentially a simple slide-rule operation, quickly accomplished when one 
mes familiar with the “feel” of the magnitude of the values involved. 
niques, which are largely personal, can abbreviate the consequent “as- 
omical values” and can also provide judgement in retaining the significant 


‘es of the computation. 
Hydr. Engr., U.S. Geol. Survey, Los Angeles, Calif. 
‘ 
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In the Appendix of the author’s paper the expression for friction head lo. 
hs, is written as 


he = £ + by, - hyo - ht 
where 


ht =m (hy, 3 hyo) 


The authors neglect to state, however, that the quantity in parentheses m 
always be positive. For an expanding reach, hy, > hyo> the quantity in pare 
theses is positive, and 


ht =m (hy, = hyo) 


Thus 
ht = f + hy, - hy, - m (hy, - bys) 
or 
Wes ag oe (hy, - hy4) (1 - m) 
If-m-= 0.5 


hy =f + (by, - hy,) (0.5) 


However, for a contracting reach, hy, < hy,, and the quantity in parenthes 
is negative. It is made positive in the following way: 


hy s=- mM (hy, = hyo) 


Then 
hg =f + hy, - hyo + m (hy, - hyp) 
- or 
hg = £ + (hy, - hy.) (1 + m) 
If m = 0.1 


hg = f + (hy, - By) (1-1) 


Expressing an approximation of a in terms of section properties, K and 
is essential to a direct method of determining the discharge of a stream 
the slope-area technique. Consider the two-section reach 1-2 in Fig. D1. | 
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Section Section 
ue 
FIG. D1. 
ym continuity 

1 1 
Qy = Ky $42 = Ko S92= Qo 
ee ae 
192 = (Ky Ky)? (S12 S22)2 


-is 


hg = £ + hy, - by. - ht 


‘a contracting reach 


ht = - m (hy, - hyo) and hp = f + (by, - by,) (1 + m) 


‘an expanding reach 


he =m. (hy, - hy») and he = f + (hy, - hy») (1 - m) 


reover 


2 2 
L QL 
he = So L= | O - Ky Ko 
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then 


Q2 L . : 
fog ee hyo) (1 + m) 


which by algebraic manipulation becomes, for a contracting reach, 


Modification of the area relationship by a values determines whethe: 
reach is hydraulically contracting or expanding; that is, whether the weigh 
velocity head is increasing or decreasing in the downstream direction. If 
quantity in brackets is positive, the reach is contracting; if negative, the re 
is expanding. 

Using data from the author’s example, in the formula for a contract 
reach, the discharge was computed as 10,700 cfs. 


PERFORMANCE OF FLOOD PREVENTION WORKS 
DURING THE 1957 FLOODS 


Discussion by Fred W. Blaisdell 


‘RED W. BLAISDELL, ! F. ASCE.—Mr. Moore has presented interesting 
rmation on the operation of the SCS flood prevention structures in Texas, 
thoma and Arkansas during the floods of April, May and June, 1957. Ac- 
ling to Mr. Moore, the performance of the structures, designed for a 50- 
torm and subjected to greater than 100-yr storms, “... exceeded expec- 
ms...” The paper is devoted primarily to the performance of the emer- 
yy spillways. This is well because it is noted that “. . . the use of vege- 
d earth spillways for flood prevention structures .. . is a material diver- 
-e from the accepted principles employed by the engineering profession in 
design.” 

egarding the principal spillway, Mr. Moore notes that it is used to dis- 
‘ge water “...ata rate designed to prevent flooding of the stream chan- 
selow.” A section through a principal spillway is shown in Fig. 3. It con- 
3 of a drop inlet, a barrel through the dam and a cantilevered outlet. The 
that the performance of the principal spillways is not mentioned is prob- 
due to the fact that they performed as designed. Thewriter has been con- 
ed with the hydraulics of the principal closed conduit spillways since 
, and naturally was extremely interested in seeing how they came through 
e severe tests. In September, 1958, he visited a number of the structures. 
visits confirmed the supposition that the structures did their intended job 
complete satisfaction. 

he practice of omitting stilling basins and simply cantilevering the pipe 
nd the toe of the fill, as shown in Fig. 3, is a practice that is also “a ma- 
il divergence from accepted principles.” It was the performance of these 
ilevered pipe outlets that the writer was particularly interested in evalu- 
r, Compared to the size of the pipe, the storage between the crests of the 
cipal spillways and the emergency spillways is large. At many of the 
ctures, this storage pool was partially drawn down and refilled several 
S as a result of successive rains and many of the principal spillways ran 
pacity for long periods of time. This gave maximum opportunity for the 
lopment of the scour holes anticipated to form at the pipe exit. 

he performance of the cantilevered outlets was much better than the 
er had expected. Typical scour holes are shown in Figs. Bl and B2. 
B1 is the outlet at Honey Creek Site No. 8E. There isa 17-in. concrete 
through the fill and the cantilevered pipe is 18-in. corrogated metal. The 


eae. |. for) — Peete ke eed ler ae ee ee 


October, 1959, by Charlie M. Moore. 
Proj. Supervisor, Agricultural Research Service, U.S. Dept. of Agric., St. Anthony 


3 Hydr. Lab., Minneapolis, Minn. 
% 
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FIG, Bl. 


FIG, B2. 
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ran for six weeks at full capacity. The 12-in. outlet pipe shown in Fig. 
for the Honey Creek Site No. 14. The depth of the pool in the scour hole 
ut 4 ft. It can be seen in Figs. Bl and B2 that the scour holes are small 
>. At both sites, as well as all other sites seen, the scour hole termi- 
close to the endof the pipe. The bent shown in Fig. 3 is intended to sup- 
he pipe if the scour hole progresses upstream. The experience of the 
structures indicates that this feature may not be necessary under the 
ions existing there. However, it is not meant to imply that bents would 
ecessary in all cases. 

2 conclusion to be drawn from the Texas experiences is that a canti- 
1d outlet may be used in place of a stilling basin. 

> writer has been looking at cantilevered outlets throughout the Midwest 
qumber of years with the expectation of finding some structures in which 


FIG. B3. 


our hole had reached a dangerous size. He recalls only perhaps two 
ures in which the scour holes looked as though they might eventually 
yp to a damaging size but which hadnot yet progressed to that point. One 
se outlets is shown in Fig. B3. This is a 42-in. cantilevered outlet built 
0. It carried a designed flow in 1952. This single flow formed the wide 
nole shown in the picture taken in 1954. Comparison with a picture 
in 1952, shows no visible increase in scour hole dimensions, possibly 
e additional capacity flows have not passed the structure. Both pic- 
show that the scour under the end of the pipe is insufficient to endanger 
ny way. In fact, the scour has not undermined the pipe as far as the bent 
to support the pipe over the scour hole. In spite of the rather wide 
10le, the outlet performance must be classed as satisfactory. Whether 
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or not future capacity flows will increase the size of the scour hole to da: 
ous proportions is of course unknown at the present time. 

The writer feels unsure of the performance of cantilevered outlets in 
of the fact that he has been unable to discover field evidence to back vu 
doubts. He feels there must be limitations to the size of cantilevered o1 
that can be used with soils of different erodibilities in order to preven 
scour hole from endangering the dam. Here isa fertile field for rese: 
Will the scour holes develop laterally, rather than in a dangerous upst 
direction? Will upstream development of the scour hole take place afte 
lateral development has reached acertain extent? How long will it take t 
velop the scour hole? What is the effect of the bed material? What will } 
ultimate size? Is this an economical self-formed stilling basin? Thes: 
questions that need answering. 

Evidence collected by the writer indicates that this type of outlet stru 
has a definite usefulness. Lack of difficulty with it is such that there : 
pressure from designers to conduct a thorough investigation into the cha 
teristics and limits of applicability of the cantilevered outlet. Howeve1 
writer feels that research on the cantilevered outlet would be of major 1 
to those who design and use this economical outlet structure. 


DISCUSSION 85 


THE VORTEX CHAMBER AS AN AUTOMATIC 
FLOW-CONTROL DEVICE? 


Discussion by Michael Amein 


ICHAEL AMEIN,! M. ASCE.—The formation of vortices in flow, in hy- 
ic structures, has been found undesirable because of three principal 
res associated with flow containing a vortex. The first is the instability 
» head-discharge relationship, resulting in an erratic hydraulic behavior. 
second feature is the reduction of the discharge capacity, resulting in in- 
ed cost for the structure. The third feature is the possibility of damage 
the erosive forces of the vortex motion. Such forces have been known to 
‘damage to inlet drop spillway structures at earthdams. The authors 
described a case in which the presence of a vortex is not accompanied 
y of the undesirable features mentioned. They have experimentally 
nstrated that the erratic behavior can be eliminated by stabilizing the 
x ina vortex chamber. In the applications envisaged by them, the ero- 
forces are either insignificant or irrelevant, and the reduction in ca- 
r would be used to advantage. 
e type of vortex described by the authors usually results from flow 
sh an orifice in the bottom of a tank, and is the closest approximation to 
»vortex. Ina free or mathematical vortex, the product of the tangential 
ity and the distance from the core isconstant. This conclusion can be 
2d by the application of Bernoulli’s theorem along a streamline, and by 
ng the radial pressure differential to the centrifugal force. In this 
3is the effects of the viscous forces are neglected. However, since ina 
luid vortex, viscous forces are present, the free vortex theory is not 
ly applicable. Furthermore, from the principles of the free vortex mo- 
s enunciated by Helmholz, the free vortex consists of the same fluid 
les. In the real fluid, the fluid particles travel in continuously shrinking 
around the core and are discharged through the orifice. Thus, the par- 
constituting the vortex are constantly being renewed. It is the pre- 
ance of the tangential velocity in the vortex chamber that gives this 
-a free configuration and behavior, so that the application of the theory 
> vortex in this case is probably correct for all practical purposes. 
2 authors have presented the results of their experiments in terms of a 
sient of discharge and adimensionless parameter called the vortex num- 
This parameter contains the prime factors in the vortex motion. The 
would favor the publication of the observed data in conjunction with the 
s. A better understanding of the physical phenomena, the precision of 


scember, 1959, by R.C. Kolf and P.B. Zielinski. — 
ydr. Engr., Anderson-Nicholas & Co., Boston, Mass. 
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measurements and the quality of experimentation could be obtained fre 
observed data. Although the presentation of data in dimensionless paran 
is a commendable procedure following the widespread acceptance of d 
sional analysis, yet, this type of presentation, unless supplemented b 
served data, freezes the data to a single interpretation. Furthermor 
data would not be accessible to support or refute other approaches that 
come to light in the future. 

The reduction in the discharge capacity of the orifice with the format 
the vortex can be said to be due to the conversion of part of the static h 
the tangential kinetic energy in the vortex. In the vortex chamber, the 
duction of flow at a tangential direction stimulates the formation of ay 
and tends to make the vortex a stable feature of the flow. The author’ 
periments indicate that the coefficient of discharge C is a unique functi 
the vortex number, from which the stability of the head discharge relatic 
is deduced. 

In the following analysis, three new symbols are used in addition 1 
ones described by the authors: b = width of either of two inlet chal 
W = velocity in the inlet channel at the entrance to the chamber; and @ = 
between the inlet channel alignment and the tangent to the chamber at the 
of entry. 

The velocity in the inlet channel is so small that the depression « 
water surface due to the velocity head is neglected. Then inflow 


The circulation is 
P68 BOWoC0s 6 ‘ 
But, from the definition of vortex number 
T=DV2¢H NG Re ea : 


Therefore, 


7B cos @ a Pee wee oe eee ew ewe 
By substituting this value of W in Eq. ls 


9 2b DE 3/2 y 
k, 7 B cos @ 


The outflow is expressed as 


o ¢ 6 8 0 © wm 


From Eqs. 5 and 6, it follows that 


7 BA cos @ 


; C _ 
He= 2b D Hidite -oonenehie: set mam 
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indicates that the head H ina given vortex chamber can determine the 
¢ number and the discharge. For the case in which b = 6 in., B = 42 in., 
in. and @ = 20°. 


C 
H = 2.1 v eiied eae Mek tats nol sids ogee. reais (8) 


en the values of vies 0.8 and Vie 2.0, the relation between C and V, as 
nined from Fig.2 is 


_ by substitution in Eq. 8, results in 


2.1H 
a i4c3inicc°: aS arte (10) 


s known, then C and Q can be readily calculated. Such a calculation is 
nted in Table 1. 

is interesting to compare the 

head required to provide a given 


\rge when there is vorticity in TABLE 1.—COMPUTATION OF 

ow to the static head required to DISCHARGE 

le the same discharge when the 

ity is eliminated. For the latter il eee es, ba ar O,; in cublone’ 
a constant C equal to 0.686 is jn feet c feet per second 


ied. The static head required to 


le a discharge Q can be deter- 0.8 0.43 0.15 

am 1.0 0.47 0.18 

1.2 0.49 0.21 

1.4 0.51 0.24 

e Q2 1.6 0.53 0.26 
2C2 A2¢ 


sults of the computations are shown in Table 2. 
s seen from Table 2, that as the head on the orifice increases, a small- 
tion of it is converted to vortex motion. Therefore, the vortex chamber 


TABLE 2 


H 
with 
vortex, 
in feet 


Be cheiis the percentage of the total static head lost in vorticity. 
ce SS —————— 


q 
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may not be useful as a flow control device, that would maintain a consta 
a slowly rising discharge under a rising head. In fact, as has been e 
enced in drop-inlet spillway structures, a sufficiently high head will c 
out the vortex motion and re-establish the normal orifice flow relations 
Many cases could exist in which the vortex chamber would perform no k 
than an orifice of smaller area. It appears that it is the range of head 
discharges within which flow takes place that determines whether the ] 
ence of a vortex in the flow can be permitted and whther the vortex ca 
made to perform a useful function. Further investigation might indicat 
even in hydraulic structures it may be more feasible, in some case 
stabilize the vortex than to make attempts towards its elimination. Resi 
on vortex motion would not only open new fields for its application, but 1 
also add to our basic understanding of the mechanics of flow. 
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(EM3), 2097(HY')°, 2098(SA4)°, 2099(EM3)°, 2100(AT3), 2101(AT3), 2102(AT3), 2103(AT3),2104(AT3), 
2105(AT3), 2106(AT3), 2107(AT3), 2108(AT3), 2109(ATS), 2110(AT3), 2111(AT3), 2112(AT3),2113(AT3), 
*414{AT3), 2115(AT3), 2116(AT3), 2117(AT3), 2118(AT3), 2119(AT3), 2120(AT3), 2121(AT3), 2122(AT3), 
2493(AT3), 2124(AT3), 2125(AT3). 


GUST: 2126(HY8), 2127(HY8), 2128(HY8), 2129(HY8), 2130(PO4), 2131(PO4), 2132(PO4), 2133 (PO4), 2134 
(SM4), 2135(SM4), 2136(SM4), 2137(SM4), 2138(HY8)°, 2139(PO4)° 2140(SM4)°. 


“TEMBER: 2141(CO2),2142(CO2), 2143(Co2), 2144(HW3), 2145(HW3), 2146(HW3), 2147(HY9), 2148(HY9), 
2149(HY9), 2150(HY9), 2151(IR3), 2152(ST7)°, 2153(IR3), 2154(IR3), 2155(IR3), 2156(IR3), 2157(IR3), 2158 
[R3), 2159(IR3), 2160(IR3), 2161(SA5), 2162(SA5), 2163(ST7), 2164(ST7), 2165(SU1), 2166(SU1), 2167(WW3), 
2 ), 2169(WW3), 2170(WW3), 2171(WW3), 2172(WW3), 2173(WW3), 2174(WW3), 2175(WW3), 2176 
‘ww3), 2177(Wws3), 2178(CO2)°, 2179(IR3), 2180(HW3)*, 2181(SA5)°, 2182(HY9)°, 2183(SUL)°, 2184 
(ww3)°, 2185(PP2)*, 2185(ST7)°, 2187(PP2), 2188(PP2). 


SER: 2189(AT4), 2190(AT4), 2191(AT4), 2192(AT4), 2193(AT4), 2194(EM4), 2195(EM4), 2196(EM4), 
(EM4), 2198(EM4), 2199(EM4), 2200(HY10), 2201(HY10), 2202(HY10), 2203(PL3), 2204(PL3), 2205 
3), 2206(PO5), 2207(PO5), 2208(PO5), 2209(PO5), 2210(SM5), 2211(SM5), 2212(SM5), 2213(SM5), 2214 
), 2215(SM5), 2216(SM5), 2217(SM5), 2218(ST8), 2219(ST8), 2220(EM4), 2221(ST8), 2222(ST8), 2223 
T8), 2224(HY10), 2225(HY10), 2226(PO5), 2227(PO5), 2228(PO5), 2229(ST8), 2230(EM4), 2231(EM4), 
(AT4)C, 2233(PL3)°, 2234(EM4)°, 2235(HY10)°, 2236(SM5)°, 2237(ST8)%, 2238(PO5)°, 2239(STS), 2240 
. 


[BER: HY11), 2242(HY11), 2243(HY11), 2244(HY11), 2245(HY11), 2246(SA6), 2247(SA6), 2248 
), so Reps tet Dr hes 2252(SA6), 2253(SA6), 2254(SAG), 2255(SA6), 2256(ST9), 2257(ST9), 
{8(ST9), 2259(ST9), 2260(HYi1), 2261(ST9)°, 2262(ST9), 2263(HY11), 2264(ST9), 2265(HY11), 2266(SA6), 
: 6), 2260(HY¥11)°, 2270(ST9). 
; eon cae Leen 2273(HW4), 2274(HW4), 2275(HW4), 2276(HW4), 2277(HW4), 2278 
), 2279(HW4), 2280(HW4), 2281(IR4), 2282(IR4), 2283(IR4), 2284(IR4), 2265(PO8), 2286(PO6), 2287 
2288(P06), 2289(PO6), 2290(PO6), 2291(PO6), 2292(SM6), 2293(SM6), 2294(SMB), 2295(SME6), 2296 
Da97(wwa), 298(WW4), 299(WWw4), 2300(WW4), 2301(WW4), 2302(WW4), 2303(WW4), 2304(HW4), 
305($T 10), 2306(CP2), 2307(CP2), 2308(ST10), 2309(CP2), 2310(HY12), 2311(HY12), 2312(PO8), 2313(P06), 
JA(ST10), 2315(HY12), 2316(H¥12), 2317(HY12), 2318(WW4), 2319(SM6), 2320(SM6), 2521(ST10), 2522 
T10), 2323(HW4)°, 2324(CP2)¢, 2925(SMé)c, 2326(WWw4)c, 2327(IRA)C, 2328(PO6)°, 2329(ST10)°, 2330 
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he 1), 2333(EM1), 2334(EM1), 2335(HY1), 2336(HY1), 2337(EM1), 2338(EM1), 
MVD, ZHOU, 2041GGai),, 2942(EM), 2840(SAL), 2944(STI), 2345(8T}, 2940(STI), 247TH), 
M1)°, 2349(HY1)¢, 2350(ST1), 2351(ST1), 2352(SA1)°, aati Da abies acveaen 

Voaseic 7(CO1), 2358(CO1), 2359(CO1), 4 : k 

), 29eA “2305601 “230082, AAOTEUL, 2368(SM1), 2369(HY¥2), 2370(SU1), 2371(HY¥2), 
, 2), 2375(PO1), 2976(HY2),,2377(CO1)°, 2378(SU1), 2379(SU1), 2380(SU1), 
, 2984(s72), 2985(8U1)°, 2366(SU1), 29871501), 2388(SU1), 2969(SM1), 


2(PO1)e. fi oe , 
2395(IR1), 2396(IR1), 9e(IR1), 2399(IR1), 2400(IR1), 2401(1R1), . 
1 Sayer an re ania Svieeh, 2410(ST3), 2411 


3 1), 406(1R1), 1), 2417(HW3), 2418{(HW1)¢, 24 awwi)?, 


' HY3), 2 
lira) BARA), Se 
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